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Abstract 

We review the X-ray spectra of the cores of clusters of galaxies. Recent high resolu¬ 
tion X-ray spectroscopic observations have demonstrated a severe deficit of emission 
at the lowest X-ray temperatures as compared to that expected from simple radia¬ 
tive cooling models. The same observations have provided compelling evidence that 
the gas in the cores is cooling below half the maximum temperature. We review 
these results, discuss physical models of cooling clusters, and describe the X-ray 
instrumentation and analysis techniques used to make these observations. We dis¬ 
cuss several viable mechanisms designed to cancel or distort the expected process 
of X-ray cluster cooling. 
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1 Introduction 


Observations show that the X-ray emission from many clnsters of galaxies 
is sharply peaked aronnd the central brightest galaxy. The inferred radiative 
cooling time of the gas in that peak, where the temperatnre drops to the 
center, is mnch shor ter than the age of the clnster, snggesting the existence of 
a cooling flow there (jFabian et al.lll994ll . X-ray spectroscopy over the past 5 yr 
shows that the temperatnre drop toward the center is limited to abont a factor 
of three. Jnst when the gas shonld be cooling most rapidly it appears not to be 
cooling at all. This is sometimes known as the cooling flow problem. Carefnl 
observations show that gently distribnted heat is reqnired over a radins of np 
to 100 kpc to balance radiative cooling in these regions. 


The issnes of cooling and heating of hot gas have broad relevance to the 
gaseous part of galaxy formation and evolution. Brightest cluster galaxies 
(BCG) are the most massive galaxies known. Calculations of the clustering 
behaviour of cold dark matter predict a power-law mass distribution for large 
galaxies whereas th e stellar mass observed has an exponential distribution 
( Benson et al The truncation of the stellar mass distribution in massive 

galaxies is likely due to the process which stops cooling flows. Simple cooling 
flows are an ingredient of semi-analytical models for galaxy formation. The 
cooling of hot gas to form stars is essential for the growth of massive galaxies 
and cannot be studied directly for isolated systems due to Galactic absorption. 
The cores of galaxy clusters offer examples which can be directly observed. 
However they do not appear to operate in any simple manner. The problem 
appears to be widespread, from the most massive clusters to the centers of 
individual elliptical galaxies. Heating and cooling problems of hot gases are 
common in astronomy, with examples ranging from the interstellar medium 
of our own Galaxy to the Solar Corona. 


The diffuse hot ionized plasma in clusters is magnetized w hich means that 
MHD processes may be important (jSchekochihin et al .11200411 . 


Here we briefly review the main X-ray properties and emission processes of 
the intracluster medium (ICM) before showing the X-ray spectra of cool cores. 
We then discuss the main solutions which have been proposed for the cooling 
flow problem. 
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2 Clusters of Galcixies 


Clusters are the most massive bound and quasi-relaxed objects in the Universe. 
They have total masses of 10^^ to above 10^® Mq. The total gas fraction is 
about 16 per cent with about 13 per cent in the hot ICM and the remaining 
3 per cent in stars in the cluster galaxies. The remaining 84 per cent of the 
mass is in dark matter. Gas densities in cluster centers range from as much as 
10“^ cm“^ in peaked clusters to 10“^ cm“^ in the non-peaked ones. This is in 
stark contrast to the mean cosmic density of baryons of about 10“® cm“^. 



Fig. 1. Chandra X-ray (left) and DSS optical (right) image of the relaxed massive 
galaxy cluster, Abell 2029. Both images are 4 arcminutes on a side. Abell 2029 is an 
extremely regular and putative cooling flow cluster. The X-ray image demonstrates 
how the intracluster medium pervades the space between the galaxies shown in 
the optical image. Figure adapted from http://www.chandra.harvard.edu/ (X-ray: 
NASA/CXC/UCI/A.Lewis et al. Optical: Pal.Obs. DSS). 

The characteristic or virial radius, i?v, of a cluster, dehned from the theory 
of structure collapse in an expanding universe as where the mean density of 
the cluster is 200 times the critical density of the Universe (i.e. 200x3iL^/87rG, 
with the Hubble constant at redshift z varying as H/Hq = ^^^(l -|- zY -|- 1 — Om), 
is typically between 1 and 3 Mpc. The gas is heated by gravitational infall to 
temperatures close to the virial temperature kT ~ GMmp/R^, which ranges 
in clusters from 1-15 keV. The total X-ray luminosities range from about 
10^3 g-i ]^g46 gj^g g-i_ Objects at lower masses and luminosities are 
groups which have from a few to tens of member galaxies as compared with 
the hundreds of galaxies in a typical cluster. 

Structure formation in the Universe proceeds in a hierarchical manner with the 
most massive objects, clusters, forming last, which means now. They continue 
to evolve by the infall of subclusters. The time since the last major merger is 
typically about 5 Gyr. About 20 per cent of clusters have had a more recent 
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merger or are undergoing one. These are not the subject of this review. 


Analytic and numerical simulations of cluster formation indicate that the total 
X-ray luminosity Lx oc in the absence of gas cooling and heating. This 
follows since the X-ray luminosity is dominated by thermal bremsstrahlung 
so Lx oc the mean gas density n oc is constant and T = ^. 


The temperature drops monotonically outward (by a factor of up to about 2). 
Observations instead show Lx oc over the temperature range 2-8 keV with 
a wide dispersion at lower temperatures and a possible flattening above. The 
simplest explanation for this result is that the gas has ha. d additional heating 
of 2-3 keV per particle ( Wn et ah 20001: Voit et ah 2003h . The effect of such 
heating is not to increase the temperature by that amount but mostly to 
expand the gas (reducing its density and thus X-ray luminosity). Such energy 
is plausibly due to energy output from active galaxies, i.e. accreting black 
holes in cluster galaxies. Alternatively, radiative cooling by rem oving the low- 
entro py gas in star formation may reproduce the relation as well (jVoit fc Brvan 


200111 . 


The gas has generally been enriched to 0.3 of the Solar value by early su¬ 
pernovae. In relaxed clusters the potential and gas peak on the BCG. The 
metallicity often rises to solar or even higher around the BCG, probably due 
to SN la. 


In relaxed. X-ray peaked, clusters the temperature prohle is often inverted in 
the inner core (i.e. R < 100 kpc) dropping inward as T oc r" with a ~ 0.3 —0.5. 
The gas density there rises as n oc r~^. 


The overall prohles of the gas density and temperature depend on the entropy 
of the gas and thus on its heating and cooling history, subject to the equation 
of hydrostatic equilibrium. 


dp dinkT) GM(< r) 

— = -pg or —--= -nprrip -r- 

dr dr r^ 


( 1 ) 


where p is the pressure, p is the mass density, n is the number density, k is 
Boltzmann’s constant, T is the temperature, G is Newton’s constant, g is the 
gravitational acceleration, M(< r) is the enclosed mass within a radius r, mp 
is the proton mass, and prup is the mean mass per particle. This equation is 
used to estimate the total mass prohle of clusters. Massive ones can act as 
gravitational lenses for background galaxies as seen in the optical band which 
provides another means to me asure mass prohles . Agreement between prohles 
determined by both methods ( Allen et al.ll2nnici1 show that hydrostatic equi¬ 
librium holds well in the main body of relaxed massive lensing clusters and 
that any non-thermal pressure there is not dominant. 
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3 Physics of the Intracluster Medium 


The intracluster medium (ICM) is plasma that is nearly fully ionized due to 
the high temperatures created by the deep dark matter gravitational potential. 
Hydrogen and helium, for example, are fully stripped of their electrons. Heav¬ 
ier elements have retained only a few of their electrons in this hot medium. 
In addition to free electrons and ions in the plasma, electromagnetic radia¬ 
tion, which is emitted mostly as X-rays, is created by quantum mechanical 
interactions in the plasma. 

The physics of the ICM can be studied as two physical phenomena: 1) the 
ionized plasma, and 2) the radiation emission processes. The ionized plasma 
is well-described by magneto-hydrodynamic theory on large enough scales. 
The radiation emission processes are governed by equations for X-ray emission 
from a collisionally-ionized plasma. We can treat these phenomena separately 
because the intracluster medium is mostly optically-thin, (i.e. the radiation 
almost completely escapes without interacting with the plasma). Later in the 
next section, we will show that this is the case. 

In §3.1, we describe how X-rays are produced in the intracluster medium. In 
section §3.2, we motivate magneto-hydrodynamic theory as the description of 
the plasma. In §3.3, we unify the two theoretical descriptions to produce what 
we refer to as the “standard cooling-flow model”. This model will serve as the 
basis for the physical description of the cores of clusters of galaxies. In later 
sections, we describe the puzzling observations that seem to agree with many 
expectations of this model, but disagree strongly with other aspects of the 
model. 


3.1 X-ray Emission from Collisional Plasmas 


Ionized plasmas produce copious amounts of X-rays. The emission of X-rays 
has two important consequences. First, it allows us to observe the intracluster 
medium by detecting those X-rays. Furthermore, since most of those X-rays do 
not interact between their emission and their detection in X-ray telescopes it 
allows us to study the intracluster medium in an unperturbed state. Through 
X-ray spectroscopy and imaging, we can measure several physical quantities, 
such as the temperature and density, at various positions in the cluster. 

The second important consequence is that the emission of these X-rays will 
tend to cool the plasma. A signihcant quantity of energy is carried away by the 
X-rays as they escape the cluster. The emission of these X-rays was thought to 
set up a non-linear process of excessive cooling in cluster cores that is broadly 
termed a “cooling flow”, a major subject of this work. In the following, we 
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describe the process of X-ray emission from collisionally-ionized plasmas and 
show how it relates to the intracluster medium. 


3.1.1 Coronal Approximation 


The emission of X-rays from ionized atoms in a plasma can be quite complex. 
Fortunately, we can make use of several approximations to simplify the emis¬ 
sion processes. The em ission of X-rays in the IC M is well-described by the 
coronal approximation (Elwert 19521: Mewel 1999ll . These approximations de¬ 
scribe optically-thin plasmas in collisional equilibrium. Collisional equilibrium 
occurs when electron collisional ionization processes are balanced precisely 
by recombination processes. The coronal approximation, as the name implies, 
was originally developed for study of the Solar Corona, but the condition also 
applies to gas in clusters of galaxies, as well as hot gas in elliptical, starburst 
galaxies, and older supernovae remnants. In this approximation, there are 
three important conditions that specify the thermodynamic state of the free 
electrons, ions, and photons in the plasma, as well as the electron distribution 
within each ion. 


The hrst approximation is that the photons are assumed to be free and do not 
interact with either the electrons or the ions after they are created. This has 
the importance consequence that photo-ionization processes (ionizing atoms 
by photons) and photo-excitation processes (raising an electron in an atom 
to an excited level) are far less frequent that electron collisional ionization 
and excitation processes. The radiation densities are low enough in clusters of 
galaxies for this condition to be met, except possibly for resonant scattering i n 
some strong emission lines with high oscillator strength (jGilfanov et al.lll987ll . 


The second approximation is that the atoms can be treated as if their electrons 
are all in their ground state rather than having a Boltzmann distribution as is 
common in LTE (local thermodynamic equilibrium) gases. This is true if there 
is a low enough electron and radiation density, so that the excitations that are 
density dependent are less frequent than the radiative decays, which are quite 
fast for X-ray energies. The radiation density is low as discussed above. This 
condition is also met for electron densities below 10^° cm“^ for even slowly 
decaying metastable states. Densities in clusters are at most 10“^ cm“^. 


The hnal approximation is that the plasma is locally relaxed to a Maxwellian 
distribution around a common electron temperature, T. The free electrons 
and ions are assumed to have obtained a common temperature. This is only 
valid if typical dynamical time-scales, such as the time it takes the plasma to 
cool, is much longer than the time scale for sharing energy between electrons 
and ions, such as the electron recombination time scale or time scale between 
Coulomb collisions. If this is true, then collional equilibrium is achieved in 














which ionizations are balanced by recombinations. This assumption is valid in 
the cores of cluster, but may break down in the outer regions of clusters where 
the density is lower and the plasma may still be in the process of ionizing. 
One might also worry that this assumption may break down in the cores of 
clusters if there is a complex multi-phase distribution of plasma temperatures 
due to thermally unstable pockets of plasma each with their own temperature. 
Generally if such a situations exists, however, it is likely that each pocket will 
locally achieve a common temperature. An electron penetrating a cloud of a 
different temperature will interact and achieve a Maxwellian distribution. It is 
only with complex and rapidly mixed interfaces that the time scale of injection 
of a new cloud is much shorter than where Ve is the electron velocity, and 
Ae is the electron mean free path, that the collisional equilibrium would be 
violated. 


3.1.2 Ionization Balance 


The fraction of atoms at a given charge state is determined by a balance of 
ionizations and recombinations. The exact balance is fixed by a coupled set of 
differential equations that relate the ionization and recombination processes 
between neighboring charge states. The equations are of the form. 


dni 

dt 


Ij^n^ni Rin^n^ R—^n^ni—x 


( 2 ) 


where rig is the electron number density, n* is the density of atoms in the Ah 
charge state, A is the ionization rate out of the Ah charge state, and Ri is the 
recombination rate out of the Ah state. We have ignored spatial diffusion of 
the ionization balance. Note, that it is customary to express the density of a 
given ion relative to the number density of hydrogen, hh- We then define a* 
as the relative abundance and the fraction of atoms in a given charge state as 
/j. Then we have, nj = where nna-i drops out of the above equations. 

Collisional equilibrium assumes that a steady-state has been achieved, so that 
the left hand side is set to zero. The equations therefore simplify to equations 
of the form. 


0 — —lift + Ri+lfi+l 


(3) 


Collisional equilibrium will eventually be achieved if a plasma remains undis¬ 
turbed for a long enough period of time (i.e., the inverse of the recombination 
rate). Several calculations have been done to determine this ionization balance 
as a function of the electron tem perature ([Jordan 19691: Arnaud fc Rothenflng 

Mazzottall998ll. One such calculation is shown 


as a mnction oi tne electron tem p< 
ll985l:lArnaud fc Ravmondl 1992l: E 
in Figure 2. 
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Fig. 2. The charge state abundance (elemental abundance times fraction ionic 
abundance) of various ions as a function of temperature. The top panel shows he¬ 
lium-like and hydrogen-like charge states of various low Z atoms. The bottom panel 
shows iron ions having the outer electron in the K, L, and M shell. The bottom 
panel indicates how the measurement of various ions in the iron series is a sensi- 
tive probe of whether plas ma at a given temperature exists. Figure uses data from 
Arnaud fc Raymond 19921. 


Once the ionization balance is determined the X-ray spectrum can be cal¬ 
culated by considering the various radiation processes. The most important 
processes are bremsstrahlung and the K and L shell transitions for the discrete 
line emission. 


3.1.3 Bremsstrahlung and other Continuum Processes 

When Hydrogen is ionized above temperatures of 2 x 10^ K, copious amounts 
of bremsstrahlung emission are produced. Bremsstrahlung radiation results 
from the accelerations of the free electrons in the Coulomb held of an ion. The 
spectrum is roughly independent of energy below the energy equal to kTg, 
where Tg is the electron temperature and k is Boltzmann’s constant. A rough 
approximation of the power per energy per volume radiated by bremsstrahlung 
is given by the equation below, 


d^P 

dVdE 




_1 E 

He Uff T 2e kT 


—— 1 
cm s 


(4) 
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where rie is the electron density, hh is the hydrogen density, E is the photon 
energy, k is Boltzmann’s constant, and T is the electron temperature in Kelvin. 
The total power radiated is therefore, 

^ 10“^'^ rie riH T2 ergs cm“^ s~^ (5) 

In addition to bremsstrahlung, bound-free emission (the capture of a free elec¬ 
tron to a bound state) and two-photon emission (which occurs most frequently 
following a collisional excitation of Hydrogen to the 2s level) are also signihcant 
source of continuum radiation. They both modify the shape of the continuum 
emission. 

3 . 1.4 Discrete Line Emission 

Discrete line emission is formed by a number of atomic processes. Such line 
emission is the most important tool for the X-ray spectroscopist. The most 
important atomic processes are collisional excitation, radiative recombination, 
dielectronic recombination, and resonant excitation. Generally, the processes 
have been incorporated in a number of publically available and well-tested 
codes that are used to study collisionally-ionized spectra. 

The strength of an emission line is determined by the excitation and recom¬ 
bination rates, which are proportional to the integral of the velocity times 
the cross-section for a particular process over a Maxwellian distribution. The 
volume emissivity for a given emission line is calculated by equations of the 
form, 

e = UeUi + C^fi + (6) 


where is the rate of inner shell ionization processes, is the sum of 
collisional excitation processes, and is the sum of recombination pro¬ 
cesses. Several well tested codes have been developed to calculate the emer¬ 
gent spectrum by assuming an ionization calculation and including a set 
of excitation and recomb inations ra tes. These include the Raymond-Smith 
f R,avmond fc Smith 1977 1. MEKA (Mewe et ah 19851 198(tI: Kaastra 19921: 
Mewe et al.ll995ll . MEKAL ( Liedahl et al.ll995ll . and APEC codes f Smith et ah 


200111 . These codes are compilations of the results of more detailed atomic 


code s, which solve the Dirac equation either by the distorted wave approxima- _ 

tion (Bar-Shalom et al.l2001 : Gu et al.l2003ll or R-matrix methods ( Herrington et al, 


I995II . The number of transitions and the accuracy of the detailed processes 


limits t he results. Extensive labo ratory work has be en applied to verify wave¬ 
lengt hs ( Brown et al .Il998[l2002ll and cross-sections (|Gu et al .119991: 1 IChen et ahl 
2OO5I1 of the transitions. 
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A number of common spectral transitions occur in most X-ray spectra. In fact, 
despite the complexity implied by the above discussion there are usually only 
a couple of dozen strong transitions that are used to determine most of the 
information that can be extracted from the X-ray spectrum. Several of the 
important emission line blends are shown in Table 1. 


3.1.5 X-ray Cooling Function and Emission Measure Distribution 

The X-ray cooling function is calculated by integrating the emission from all 
processes and weighting by the energy of the photons. 

7 da 

A{T,Z,)= J dE E—{E,T,Z,) (7) 

0 


where ^ is the energy dependent line power (or continuum power). The cool¬ 
ing function relates the total amount of energy emitted per volume for a 
given amount of plasma with a given temperatur e and emissivity. The cool¬ 
ing function has been comp iled in various tables ( Bohringer fc Hensler 19921: 


Sutherland fc Dopitalh 99311 . 



Fig. 3. The radiative cooling function is shown for solar abundances (top curve), 
one-third solar abundances (middle), and pure Hydrogen and Helium (bottom). The 
X-ray region from 10® to 10® K is dominated by bremsstrahlung at high tempera¬ 
tures as well as significant contribution from line emission at lower temperatures. 
Below 10® K in the UV temperature range, the cooling function rises significantly. 
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Table 1 

Important line blends in Cluster X-ray Spectra 


Ion“ 

Wavelengths 

Energies 

Temperature 


A 

keV 

keV^ 

Fe XXVI 

1.8 

6.97 

> 3.0 

Fe XXV 

1.9,1.9 

6.70,6.63 

1.0 ^ 8.0 

Fe XXIV 

10.6,11.2 

1.17,1.11 

0.9 ^ 4.0 

Fe XXIII 

11.0,11.4,12.2 

1.13,1.09,1.02 0.8 ^ 2.0 

Fe XXII 

11.8,12.2 

1.05,1.02 

0.6 ^ 1.5 

Fe XXI 

12.2,12.8 

1.02,0.97 

0.5 ^ 1.0 

Fe XX 

12.8,13.5 

0.97,0.92 

0.4 ^ 1.0 

Fe XIX 

13.5,12.8 

0.92,0.97 

0.3 ^ 0.9 

Fe XVIII 

14.2,16.0 

0.87,0.77 

00 

o 

T 

CO 

o 

Fe XVII 

15.0,17.1 

0.83,0.73 

0.2 ^ 0.6 


15.3,16.8 

0.81,0.73 


S XXVI 

4.7 

2.62 

> 1.0 

S XXV 

5.1,5.0 

2.43,2.46 

0.3 ^ 1.0 

Si XIV 

6.2 

2.00 

> 1.0 

Si XIII 

6.6,6.7 

1.87,1.84 

0.2 ^ 1.0 

Mg XII 

8.4 

1.47 

> 0.7 

Mg XI 

9.2,9.3 

1.35,1.33 

0.1 ^ 0.6 

Ne X 

12.2 

1.02 

> 0.4 

Ne IX 

13.5,13.7 

0.92,0.90 

0.1 ^ 0.3 

0 VIII 

19.0,16.0 

0.64 

> 0.2 

0 VII 

21.6,22.0 

0.57,0.56 

0.1 ^ 0.2 

N VII 

24.8 

0.50 

> 0.1 

C VI 

33.7 

0.37 

> 0.1 


“ Note this line transition list is somewhat crude since we have chosen to tabulate 
line blends rather than actual transitions, but matches well with the quality 
of the observations. 

^ The temperature range is calculated to roughly show the temperature range 
where the emissivity of a given ion blend is within an order of magnitude of 
its peak emissivity. 
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The relative distribution of plasma at a set of temperatures is often expressed 
in terms of an emission measure distribution. The differential emission mea¬ 
sure, is dehned by 


de 

dE 


OO 


0 


dEM da 
dT dE 


{E, T) dT . 


( 8 ) 


where ^ is the energy-dependent emissivity and the integral is over all temper¬ 
atures. It is also convenient to express the distribution of plasma temperatures 
in terms of the differential luminosity, which is dehned by 


dL 

dT 


dEM 

dT 


A(T) . 


(9) 


3.1.6 Cooling Time 

The cooling time of an optically-thin plasma is the gas enthalpy divided by 
the energy lost per unit volume of the plasma. The gas enthalpy is |n/cT and 
the energy lost per volume is the electron density squared times the cooling 
function. The cooling time can then be written as, 

-nkT 

tcooi = ~ tnTsAzlsnZl ( 10 ) 

where tn is the age of the universe (13.7 Gyr), Tg is the temperature in units 
of 10® K, A _23 is the cooling function in units of 10“^® ergs cm® and n _2 
is the density in units of 10“^ particles cm“®. We used the gas enthalpy per 
volume, |n/cT instead of the thermal energy per volume |n/cT since the plasma 
is compressed as it cools which therefore effectively raises its heat capacity by 
a factor of |. Therefore, X-ray plasma with gas density above 10“^ cm“® has 
had sufficient time to cool. In the cores of cooling clusters the cooling time 
approaches cooling times below 5 x 10® yr. If the gas was nndisturbed, it would 
have a chance to cool several times. Note that as gas cools at constant pressure 
(due to the weight of overlying gas) then the rise in density as the temperatnre 
drops means that tcooi becomes shorter and shorter. 


3.1.7 Optical Depth, Resonance Scattering, and Opacity 

The optical depth for photons of a given wavelength, r, is the product of the 
column density of a particular ion, W and the cross-section of a particular 
process. The cross-section is a function of energy. The column density is the 
line integral of the ion density, Uj, which is a function of the spatial position 
along the line of sight. 
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The intracluster medium is, for the most part, optically-thin to its own radi¬ 
ation (i.e. photons escape once they are emitted without interacting again). 
Photons with energies close to certain resonance transitions, however, can 
scatter several times before leaving a cluster. The optical depth to a reso¬ 
nance transition is given by 




r = 


dlriifi 


Tie exp 


rUpC 


yiinr? 


( 11 ) 


where n* is the ion density, fi is the oscillator strength of the transition, rrie 
is the electron mass. A* is the wavelength of the transition, c is the speed of 
light, e is the electron charge, and a* is the line width given by. 



where A is the atomic mass number, p is the mean mass per particle, M is the 
Mach number of the turbulence or gas motions in the plasma, k is Boltzmann’s 
constant, and T is the temperature. The above expression includes both the 
thermal broadening of the line as well as turbulent broadening. 

If a signihcant quantity of lowly ionized matter exists along the line of sight. 
X-rays can be absorbed and re-emitted at lowly ionized longer wavelengths. 
This situation occurs frequently for absorption from neutral gas in the Milky 
Way Galaxy, but could also occur from gas trapped in the cluster potential. 
This subject is discussed in detail later, but we note that understanding the 
role of absorption often has a significant effect on the interpretation of the 
soft X-ray spectrum. This is particularly true at low resolution. Helium K 
shell absorption at low energies and Oxygen K shell absorption at 23.5 A are 
the largest contributors to the opacity from a neutral absorber and produce 
absorption edge features in the spectrum. 


3.2 Magneto-hydrodynamics 


The plasma of the intracluster medium can be described by magneto-hydrodynamics 
on large scales. Several assumptions underlie this statement. First, the ICM 
is optically-thin so we do not have to include radiative forces in the theory. 
Second, we assume the ICM is non-relativistic. This is true since the sound 
speed is at most a few thousand kilometers per second. Third, we assume the 
ICM is nearly fully-ionized, which is true since Hydrogen consitutes most of 
the plasma. Finally, the plasma can be treated as a fluid using continuous 
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fields if the plasma parameter is small so that collective processes dominate. 
The plasma parameter is defined by 


9 = 


1 

n\% 


8 X 10 ^n2T 2 


(13) 


where A/) is the Debye length, n is the density in units of cm and T is the 
temperature in K. The ICM plasma parameter is in fact extremely small. 


Magnetohydrodynamics (MHD) in this context is described by the following 
variables: the fluid velocity, v, the temperature, T, the density, p, and the 
magnetic field, B, and the dark matter density, pdm- These quantities are 
related to one another by the full set of MHD equations, which assume local 
thermodynamic equilibrium of these quantities. The transport properties of 
the fluid can be written in terms of the viscosity, p, conductivity, K, and 
resistivity, A. 


We also adopt the following notation. T is the energy per particle; whereas Tx 
is the temperature in Kelvin. They are related to one another by the relation, 
T = kTx/{pmp). k is Boltzmann’s constant, the mass of the proton is nip, 
and the mean mass per particle is pnip. Additionally, the cooling luminosity 
is usually defined in terms of the electron, Ue and hydrogen, uh, number 
densities, which are related to the fluid mass density by, Ue = nn/l.lQ = 

p/{pmp). 


Below are the 5 MHD equations. 
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The first equation, the mass conservation equation, is a statement that the 
total mass of the fluid is constant. The second equation, the Navier-Stoker 
equation, enforces momentum conservation. This equation relates the momen¬ 
tum of the fluid (left hand side) to the gravitational compression (hrst term 
on the right hand side), the thermal and magnetic pressure gradients (sec¬ 
ond and third term), tangetial magnetic transport (fourth term), and viscous 
forces (last terms). 


The evolution of the magnetic field follows the induction equation, the third 
equation. The hrst term on the right hand side generates the magnetic held 
due to plasma motion, and the hnal term dissipates the magnetic hux due to 
magnetic reconnection. 


The energy equation, the fourth equation, expresses the balance between heat¬ 
ing and cooling terms. The left hand side describes the energy content of the 
plasma (hrst and second term) as well as its compression (third term). The 
right hand side contains a conduction term (hrst term), a magnetic dissipa¬ 
tion term (second term), viscous heating terms (third and fourth term), and 
the energy lost due to radiative cooling. In addition, this equation could in¬ 
clude interactions between this plasma and other matter, such as dust, cosmic 
rays, or dark matter. Finally, the gravitational held is set by the hfth equation 
where both the dark matter and plasma contribute to the gravitational held. 
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The transport coefficients for the intracluster medium is the subject of much 
theoretical work. The values of the se coefficients f or a i onized plasma with no 
magnetic field was worked out by Snitzer et ah ( 1962ll using kinetic theory. 
These values are given by: 



A ta ngled magnetic held due to MHD turbulence, e.g. iGoldreich fc Sridhar 


I995L could modify these c oefficient considerably and there has been consider¬ 


able debate in this subjec t (iTaoll 99.4 1 IChandran fc Cowlevlll994 l iNaravan fc Medvedev! 


l2nnit iMaron et al.ll2nn4ll . We will return to this subject in later sections. 


3.3 Cooling flows 


The radiative coo ling time in the cores of at least two thirds of low redshift 

299211 


and moderate redshift ( Bauer et af 


fjPeres et ahl ll 

than 19 Gyr and for one third it is less than about 3 Gyr. The energy loss 
is directly due to the observed X-ray emission with no major bolometric cor¬ 
rection. If there is no heating to compensate the cooling then a cooling how 
occurs in these regions. In order to understand what the ’cooling-how prob¬ 
lem’ is, why heating is required, how a ’residual how’ might operate and what 
happens when heating is not operating, we now briehy examine cooling hows. 


3.3.1 Single-phase flows 

The radiative cooling time fcooi at tens of kpc radius in a cluster always exceeds 
the gravitational dynamical time so cooling leads to a slow, subsonic inhow 
there. The how causes the density to rise and so maintain the pressure, which 
is determined by the weight of the overlying gas. 

We can simplify the MHD equations signihcantly for the simple case of an 
unmagnetized single-phase subsonic how. We also ignore any terms with con¬ 
duction, viscosity, and resistivity. If we simplify the previous equations and 
rewrite the LHS of the energy equation based on the mass equation. 
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§^ + V^(pv)=0 

dvi d<^ d{pT) 

p—+ p(v.V).. = p— 


(22) 

(23) 


d 


dt 


1 f d 


— + v- v-r — T^ + v- Vp 


p \dt 
V^<I> = A-kG{p + pdm)- 


= —nenH^{T, Z) + interactions(24) 

(25) 


We now assume a spherical geometry and assume that the system is in a 
steady-state such that the partial time derivatives can be ignored, and assume 
that the flow is subsonic, then terms of order can be ignored. We combine 
the second equation with the third, and define M in the first equation, and 
obtain. 


M = pv = constant 


dA> _d [pT) 
^ dr dr 


pv— ( -T — <f> = —nennMT, Z) + interactions 

dr \2 


d /5 


1 d f 2^^ 


dr V dr 


^ = 47rG(p -h pdm)- 


(26) 

(27) 

(28) 
(29) 


When the gravitational held is important, then the temper ature approaches a 
critical soluti on which ’follows the gravitational potential’ (jFabian et ahlll984 ; 
Nnlsen 198fill as can be seen from the hydrostatic equation for a power law 


solution. Essentially the flow s ettles down to a tem perature profile close to the 
local virial one. In an NEW (jNavarro et ahl 1199711 potential where the inner 
power-law part has dark matter mass density varying as p oc r~^ this means 
T (X r. The temperature then flattens to T ~ 1 keV within the inner 10 kpc 
where the gravitational potential of the central galaxy, assume isothermal, 
dominates. The temperature Anally collapses at the center, before which the 
flow may go supersonic (the inertial velocity term is then needed in the above 
momentum equation). Over most of the region r/v ~ tcooh as seen in the 
energy equation, which varies as /n for bremsstrahlung. Using the mass 
flow equation to substitute for n, n oc r“4 in the NFW case which leads to a 
steeper surface brightness profile (the emissivity is proportional to n^) than 
observed. 


19 



















Even with the King potential commonly used before the 1990s, it was realized 
that the central sur face brightness was too steep to match the observed prohles 
fjFabian et al.lll984ll . Interpreted as a cooling flow, the data indicated that the 
mass flow rate increased with radius (roughly as M oc r out to the cooling 
radius Tc where tcooi ~ 10^° yr). Such a situation requires that matter is cooling 
out over a range of radii, which was explained by gas with a range of densities 
and thus cooling times occurring at a given radius. Therefore, more generic 
models than a simplihed spherical single phase flow had to be considered. 


3.3.2 Thermal Instability and Multiphase flows 


The above discussion of single phase flows says little about what happens on 
small scales. In particular, when cooling begins and what size perturbation 
leads to the largest growth rate and whether the multiphase flows could de¬ 
velop like the data seemed to indicate are an open question. There has also 
considerable effort to understand whether the reservoir of heat in the outer 
regions of clusters can be transferred to the center, which could effectively 
stabilize any initially thermally unstable parcels of gas. 


Field f 1965ll originally discussed the origin of thermal instability due to the 


emission of radiation. He found that for all X-ray temperatures the gas is 
thermally unstable and the growth rate is fastest on the smallest scales. He 
further found that the small scale perturbations are damped by conduction 
so that the gro wth rates will be fas t er on somewhat larger scales. A num¬ 
ber of authors (iM alagoli et ah 1987: White fc SarazinI 1987a : BalbusI 19881: 
Loewenstein Il989l: iBalbns fc S okeii Il9 89ll s tudied thermal instability in the 


context of gravitational held. iBalbna ^ 199 noted that some of the thermal 
instability arguments are inapplicable in a cluster gravitational potential. It is 
possible that overdense parcels of plasma can com e to equilibrium at a lower 


adiabat deep in the potential ( Cowie et ah 1983) • Kim fc Naravan (20033) 


argue that the radial modes are unstable even in the presence of conduction. 


NulsenI (1198011 and iThomas et ahl (1198711 . however, considered the development 
of perturbations into cooling flows and discussed the possible importance of 
magnetic helds in p inning parcels of plasma to the general hydrodynamic flow. 


LoewensteinI ()1990ll discussed the importance of the magnetic held in altering 


the instability c o nditio ns by effectively elimina ting the buoya ncy problem by 
Balbus fc Soker ( 1989ll with magnetic stresses. iBalbusI (1991) later conhrmed 


these instabilities but stressed the importance of inefficient conduction for 
these conditions. However a particular, and not explained, spect rum of density 


perturbations is required to o btain the inferred relation M oc r (jThomas et al 


1987 : Tribble 1989ll . Binnev (2004) argues that multiphase hows do not occur 


in real clusters. 
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3.3.3 Standard Multiphase Cooling Flow Model 

The standard multiphase model simplifies the physics of the flow by assuming 
that radiative cooling dominates the flow and then looking at the relative 
amount of radiation emitted at each temperature over the cooling flow volume. 
This can be seen as a starting point to testing the idea that cooling flows are 
operating in the cores of clusters. If we restore Equation 24 to the full time- 
derivative form and then integrate over volume and neglect any heating or 
additional cooling we obtain, 


We define a mass loss rate per time, M, and a differential X-ray luminosity, 
dLx, according to 


dV 


dLx = UeUndiiT, Z)dV. 


(31) 


Then the Equation 30 simplifies to. 



dLx. 


(32) 


This can be rewritten as. 


dLx 

It 



1 dp \ 
~p'dT) ■ 


(33) 


In general, ^ will be set by the local gravitational field and magnetic pres¬ 
sure. If the gravitational field is relatively smooth as can be expected in dark 
matter haloes, then the pressure will remain nearly constant in small regions 
that will begin cooling. Whether this occurs depends on the development of 
thermal instability as discussed in the previous section. If ^ term is small, 
this expression simplifies to 


dLx 5 Mk 
dTx 2 pmp 


(34) 


which is known as the standard isobaric cooling flow model. If the density 
is constant which would result from high magnetic pressure, it alternatively 
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reduces to 


dL^ 3 Mk 

dTx 2 finip 


(35) 


which is for isochoric cooling. An adiabatic equation would result if cooling 
were ineffective and the gravitational potential was strong and the left hand 
side would be zero. In general, ^ could be quite complicated, but note that 
it is fairly well established observationally that over the whole cooling flow 
volnme the temperatnre derivative, ^ is positive and the pressnre derivative 
is negative, Therefore, the last term is likely to be positive, and the relative 
amonnt of X-ray luminosity would be emitted according to somewhere in 
between the previous two equations if X-ray radiative cooling dominated the 
energy release. 

Note that Equation 33 is just the first law of thermodynamics {dQ = dU + 
pdV = ^NkdT+d{pV)-Vdp = \NkdT-Vdp = ^NkdT-Nf) differentiated 
with respect to time {,\Nk^ — ^^kT = with X-ray cooling as the only 
heat loss term. 

Equation 34 is also frequently expressed in terms of the emission measure, 

dEM b Mk I 

= o-T 

diK 2 pmp A 


This then can be used in conjnnction with the atomic physics necessary to 
produce an X-ray spectrnm as shown in Figure 4, which can be compared with 
the data. X-ray spectroscopy has demonstrated that this model is incomplete 
and therefore we have most likely neglected additional heating or possibly 
cooling in our derivation of this model. The following section describes the 
instrnmentation and analysis necessary to test it. 
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Isobaric Multiphase Cooling-Flow Model 



Fig. 4. The standard isobaric cooling-flow model. The model is produced by sum¬ 
ming collisionally-ionized X-ray spectra within a temperature range such that the 
relative amount of luminosity per temperature interval is a constant. This model 
predicts relatively prominent iron L shell transitions between 10 and 17 A that arise 
from a wide range of temperatures. In particular, Fe XVII which is emitted between 
500 eV and 1 keV has strong emission line blends at 15 and 17 A. 
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4 X-ray Instrumentation and Observational Techniques 


A number of X-ray instruments have been launched above the atmosphere, 
which is opaque to X-rays, to study X-ray sources in the last 40 years. X-ray 
telescopes typically reflect X-rays at grazing incidence using mirros coated 
with a high-Z material. CCDs, proportional counters, and microchannel plates 
have been used to record the positions of X-ray photons and make low resolu¬ 
tion energy measurements. Crystals, reflection and transmission gratings have 
been used to disperse X-rays and produce high resolution spectra. 

Below, we first discuss some of the important characteristics of X-ray tele¬ 
scopes and the relevant instrumentation that is used to study clusters of 
galaxies in X-rays. Then, we discuss the analysis techniques that are used 
to interpret the data collected from these satellites. 


4-1 X-ray Telescopes and their Relevance to Clusters 


X-ray photons from astrophysical sources have been detected by a number of 
X-ray telescopes. For observations of clusters of galaxies, several instrument 
characteristics have been important for gaining a complete description of the 
X-ray emission. 

Spectral Resolution: First, the spectral resolution or resolving power (-^ 
or -^) affects the ability to extract useful information from the spectrum. For 
typical temperatures in clusters of galaxies, the Fe L complex most strongly 
constrains the distribution of temperatures. Merely detecting or not detecting 
an emission line from a given charge state from a given ion, places narrow 
constraints on the distributions of temperatures in the plasma. The required 
resolution to resolve lines from individual Fe L ions is about 100 at 1 keV. Most 
X-ray detectors, which historically have been either proportional counters or 
solid-state devices, have been unable to achieve this resolution and therefore 
the Fe L complex appears as an unresolved bump in the spectrum. Propor¬ 
tional counters work by the incoming X-ray being photoelectrically absorbed 
an inert gas atom followed by the measurement of the electrical discharge in¬ 
duced by the cascade Auger electrons and the further release of fluorescent 
X-rays. They typically achieve resolution near 20%. CCD devices work by the 
X-ray creating a photo-electron in the silicon, which creates a electron-hole 
pair cascade through ionization. The electrons drift to a set of surface con¬ 
tacts due to an electric held and the number of electrons is used to measure 
the energy of the incident X-ray. These typically have resolving powers near 
10 or 20 at 1 keV. These devices can, however, measure an average temper¬ 
ature accurately from the shape of the continuum, which is sensitive to the 
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bremsstrahlung emission, or from the centroid of the Fe L complex. These 
instruments, however, have difficulty contraining the distribution of tempera¬ 
tures composing a single spectrum. 

Dispersive instruments, which either use a grating or crystal, have achieved 
high spectral resolution in the Fe L band. Crystals disperse X-rays accord¬ 
ing to the Bragg condition in which constructive interference is set up by the 
X-ray wavelength being equal to a multiple of the projected crystal spacing. 
Reflection and transmission gratings disperse X-rays by a constructive peri¬ 
odic surface where the X-ray wave will interfere constructively for a given 
angle of incidence. Both the Reflection Grating Spectrometers on XMM- 
Newton, and the Focal Plane Crystal Spectrometer on the Einstein Observa¬ 
tory produced high resolution X-ray spectra of clusters. Future missions may 
include non-dispersive microcalorimeters, which use the temperature change 
of a cryogenically-cooled absorber due to the photo-electron to measure the en¬ 
ergy of the X-ray. The Astro-El and Suzaku missions included a microcalorime¬ 
ter, but unfortunately the rocket failed for Astro-El and the microcalorimeter 
did not operate long enough to observe cosmic X-rays for Suzaku. A mi¬ 
crocalorimeter can achieve sufficiently high spectral resolution in the Ee K 
band such that few hundred km/s velocities of the intracluster medium could 
be measured. 

Effective Area and Exposure Time: The effective collecting area and the 
exposure time of the observation determine the number of photons collected. 
The number of photons enters into the effectiveness of the observation in 
a number of ways. Typically, tens of thousands of photons are collected in 
a given observation. Hundred of photons are necessary to make any detailed 
surface brightness image. Thousands of photons are necessary to make detailed 
spectroscopic measurements of the temperature. Tens of thousands of photons 
are necessary to make detailed measurements of elemental abundances as well 
as the construction of the differential emission measure. 

Field of View: Field of view determines the fraction of photons from a cluster 
that are detected. Often, the held of view has been smaller than the sizes of 
typical nearby clusters and therefore are an important factor in comparing 
results from different instruments. 

Spectral Bandpass: The spectral bandpass of most X-ray instruments has 
been quite large and for the most part encompassed most of the X-ray spectral 
band. Most of the photons from clusters of galaxies are emitted below 2 keV, 
and most instruments have had high efficiency at these energies. 

Angular Resolution: Angular resolution is an important factor for spatially- 
resolved spectroscopy. XMM-Newton (with FWHM of 6”) and Chandra (with 
FWHM of 0.5”), for example, have been able to perform spatially-resolved 
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photometry on kiloparsec scales for nearby clusters. 


4-2 Analysis Techniques 


Analysis techniques play a vital role in the interpretation of X-ray spectra from 
clusters of galaxies. The data analysis is usually quite complex compared to the 
techniques that can be applied to data from unresolved sources. In addition, 
the instrument response functions are actually quite complex, and therefore 
the application of these functions is usually problematic and results often are 
somewhat model dependent. 

X-ray photons that are recorded in detectors after being reflected and dis¬ 
persed by optics have three quantities that are measured: two detector coor¬ 
dinates, X and y and one energy measurement, p. These are indirectly related 
to the position of the photon on the sky, 0 and 0, and the photon’s intrinsic 
energy, e. For dispersive spectrometers the relationship between these three 
variables and their detector counterparts are very indirect and the full convo¬ 
lution has to be considered. For non-dispersive spectrometers, the relationship 
is more direct so approximations can be used. For example, it is customary to 
assume that (x, y) ~ (0,0) when extracting a spectrum to perform spatially- 
resolved spectro-photometry or to assume that p ~ e to construct an image. 

The detection probability, D, for a photon emitted from solid angle position, 
Q, and energy, e with measured values of detector coordinates (x, y) and CCD 
pulseheight p is given by 

Dix, y,p) = J de dVt R{x, y,p\ e,n) ^ (37) 


where 


d?F 


dQ.dE spatially-varying spectral source model and R is the instru¬ 

ment response function. All existing data analysis methods either solve this 
integral by, 1) making approximations abou t the response function, 2) use a 
Monte Carlo approach ( Peterson et al in®, or 3) perform an integration over 
the solid angle and just study the spectra 4) or perform an integration over 
the spectra and just study the image. One can easily see some of the reasons 


for confusing results on spectra from clusters of galaxies. It is often somewhat 
ill-dehned by what one analysis would mean when compared to an analysis 
from a different instrument without specifying the full source function, 


Note that most analyses construct a response matrix and an ancillar y response 

or SPEX 


hie, which are used in spectral codes like XSPEC (jArnandl 
( Kaastra IQQfill . These codes multiply a source spectrum on a grid by the 
response matrix, which relates the input spectrum to the model observed 
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spectrum. The model observed spectrum can be compared to the real data 
by statistical tests, typically a binned calculation. The response matrix 
{Rm) /ancillary response (A) hie approach is a simplihcation of the response 
function above by assuming the following separation, 

Rix,y,p I e,0,^/) ~ RMip \ e)A{e) (38) 


where A is constructed separately for each spectrum depending on the spatial 
position [x, y) of the photons being studied. This approach may or may not 
be appropriate for a particular analysis. 

The source function is in turn related to the function in Equation 7 in §3.1.5 
differentiated with respect to spatial coordinates. 

TT = 7 (39) 

(ladE J dE dndT ’ 

0 


where is the differential emission measure per solid angle. This function 
is the true observable if the plasma is assumed to be in collisional equilibrium. 
Note that if a given patch of the sky is assumed to be isothermal this will 
directly give a measurement of the density, since the volume element can be 
related to the angular coordinate hi for an assumed source distance, assuming 
some geometry to obtain the depth. The temperature structure in clusters 
is never completely isothermal, so a variety of temperatures will be sampled 
along the line of sight. 

We will not review all work in dealing with these data analysis problems, but 
refer to the individual work in the following section for more details since it 
can vary considerably between authors. 
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5 X-ray Spectra of Cooling Clusters 


A wealth of information has been collected by a number of X-ray satellites 
on cooling core clusters. Recently, there have been some strict tests to var¬ 
ious aspects of the cooling-flow scenario. We review many of the previous 
spectroscopic and imaging results that led to the formation of the cooling- 
flow model for cluster thermal evolution. We then discuss how high resolution 
X-ray spectroscopic results have questioned a number of assumptions in the 
interpretation of the previous work. These observations have largely galva¬ 
nized support around the idea that the thermodynamic evolution of cluster of 
galaxies is considerably different than was previously thought. We then review 
some recent X-ray imaging and spectro-photometric work that has left some 
clues to the resolution of the cooling-flow problem. Finally in this section, we 
deflne the structure of the cooling-flow problem in clusters of galaxies. 


5.1 Early Work on Imaging Observations 


The cooling time of the ICM is less than the Hubb le time at the centers 
of many clustes as demonstrated by Lea et al. (IQT^) with the Uhuru mis¬ 
sion. Imaging observations established the existence of sharp surface bright¬ 
ness peaks in some cluster of galaxies. These surface brightness peaks and 
short cooling times were interpreted as strong evidence for the existence of 
cooling flows. The existe nce of sharp surface brightness peaks w ere found by 


the C opernicus satellite (iFa bian et al' 1974 : Mhchel^^t,^ 


sions ( Gorenstein et al. 1977ll . and SAS-3 ( Helmken et al 


197 ^ 


1975 


rocket mis- 
Later mis¬ 
sions established the widespread nature of the phenomenon by establishing 
that approximately 50% of all cl usters contained sharp surface brig htness pro¬ 
files. The Ein stein Observatory (Stewart et al. 1984 : Arnandl 19981) , EXO SAT 
Observatory ( Lahav et ah 198^: Ec ge fc Stewaril 199l[ Edge et al. 199^, and 
ROSAT Observatories (|Allen et ah 2nniall confirmed this picture with obser¬ 
vations of hundreds of nearby clusters of galaxies. 


Considerable work wa s done to establish the implied mass deposition rates 
from the cooling flows (j, Tones fc Formanlll98411 . The mass deposition rates, M, 
were estimated to be between 1-10 solar mass per year for a giant elliptical 
galaxy and 100-1000 solar masses per year for the largest clusters of galaxies. 
A first order estimate of the mass deposition rate is simply made by taking 
the total gas mass calculated from the density profile and dividing that by 
the cooling time at the edge of the region. Alternatively, this can be written 
as the luminosity of the cooling flow volume divided by the temperature (in 
units of kT) times the mean mass per particle. More detailed calculations take 
into account the work done on the gas by its compression in the gravitational 


















































poten tial ( Thomas et al. 1987; White fc Sarazin 1987a 3c) 19881: Allen et ah 


I 993 I 1 . Generally, this effect is thought to reduce the overall level of mass cool¬ 


ing by at mos t a factor of tw o from the amount predicted without gravitational 
compression ( Arnand 11 


5.2 Early Work on Low Resolution Spectroscopy 


The modeling of the X-ray spectra of cooling clusters has become considerably 
more mature in 30 years. In addition, some of the early work with cooling- 
flow models has now been discounted. This work exemplifies the difficulty in 
applying models at low spectral resolution. We briefly discuss some of the 
major changes to modeling the X-ray spectrum of cooling clusters. 

The Areil-V satellite established that the emission from clusters of galaxies 
was thermal (i.e. likely in collisional equilibrium) by detecting the emission 
lines fro m the Fe Ka blend (Fe XXV Hea, Fe XXVI Lya) from the Perseus 
cluster ( Mitchell et al. 197 (tI 1 . Previously, there were arguments for the origin 
of the X-ray emission from clusters as either inverse compton scattering of 
microwave background photons off relativistic electrons or emission from a 
collection of unresolved X-ray binaries. The detection of thermal emission, 
however, established that the hot gas was trapped in the dark matter potential 
and emitted at a temperature close to what was expected from the velocity 
distribution of member galaxies. 


After the thermal model was established, hundreds, if not thousands, of clus¬ 
ters had their “mean” temperature measured with solid state detectors or 
proportional counters most frequently on the Einstein, EXOSAT, Beppo- 
SAX, and ASCA observatories. This work established relationships between 
many observables in clusters. In particular, relationships between tempera¬ 
ture and luminosity (IPavid et al. i'9931: White et al. 1997 : Markevitchl 19981: 


I 999 I: Horner et al.l I 


Wu et ah.._. . . 

20011. and entropy Ie 


199£ 


mass 


(|Nevalianen et al .1 12 OOOI: IXu et al 


oumau et al.lllOo'^ indicate that mean cluster gas prop¬ 
erties roughly follow from what is expected from self-similar gravitational col¬ 
lapse. There may be some deviation from these relations, more preferentially 
in lower mass clusters indicating the importance of non-gravitational processes 
in cluster formation and evolution. 


It was dis covered quite early that many clusters do not have an isothermal 
structure ( Ulmer fc Jernigan 1978ll and have softer X-ray spectra in the core. 
This implies lower mean temperatures. A number of more complicated spectral 
models were later used to represent the X-ray emission from cooling-flows with 
EXOSAT, Einstein, Beppo-SAX, and ASCA observations. Two temperature 
models were used by a number of authors to represent the distribution of tem- 
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Fig. 5. The X-ray spectrum of the central region of the Virgo cluster of galaxies 
obtained with the ASCA observatory. This spectrum is one of the highest quality 
spectra obtained prior to the launch of the RGS instruments on XMM-Newton. The 
bumps in the spectrum are due to unresolved blends of emission lines. In particular, 
the bump near 1 keV is due to a forest of Fe L shell transitions, which are quite 
sensitive to the temperature of the plasma. 


peratures in the cores of clusters (Matsumoto 1996 : Ikebe et et al 


19981: iFukazawa et al.lll998L l2000l: iMakishimalbnm h . An alt ernative inte r preta- 


tion of the cold gas in the cores of cluster was developed bv iMakishimal ()2nnih 
where the dark matter halo of the cD galaxy is responsible for an unmixed 
ICM component. 


Others used the cooling flow model discussed in §3.3 modified by an absorber. 
The absorb er was necessary to reduce the soft X-ray emission at low energies 
observed bv lWhite et al.l ()l991ll . In retrospect it is clear that these observations 
were demonstrating a lack of emission at low energies relative to the cooling 
flow model, although this was not explicitly stated. In particular, the relatively 
low spectral resolution made it difficult to distinguish between models. 


The a bsorber in these models was modelled typically as a single absorbing 


scree n (.lohnstone et al.| l992 ; Fabian et al.ll 994 ; Rnote X Fabian 


scree n iLyynna_eT_^Lili^uan^ 
Tiii lBuote et al.ll999ll2000al:lAllen 


2nmah . The screen was placed 


between the cluster and the observer, but the assumption was that this ab¬ 
sorption was representing cold gas in the center of the cluster. The absorption 


typically had column densities near lO^^cm which would have implied large 
quantities of cold molecular gas if it was distributed uniformly. Initial attempts 
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to find this quantity of gas were unsuccessful (lO’Dea et alJll 


alth ough 


considerable cold materi al has since been found by Edge et al.l ()2nnill and 
Salome fc Comb^ ( 2004 1. Another model used was a cooling flow model with 


a single absorption edge. The absorption edge had an energy near 0.7 keV, 
which could have been fro m ionize d Oxygen possibly from warm (10® K) ma¬ 
terial ( Bnote et ah 2000bll or dust ()Allen et ah 2001bll . Later, high resolution 
Reflection Grating Spectrometer (RGS) observations demonstrated that the 
single edge model mimicked the effect of lack of emission from colder Fe ions. 
It is interesting to note, in particular, that the htted energy of the single ab¬ 
sorption edge was slightly lower than the 3s-2p transitions of Fe XVII, and 
would therefore optimally appear to absorb cooler emission from the standard 
cooling flow model. 


Many observations generally established the cooling time is a fraction of a 
Gyr in relaxed clusters. They also clearly indicated that the temperature was 
lower where the central surface brightness was highest. These two fact are now 
shown clearly with recent data from the Ghandra observatory in Figure 6 and 

7. 


5.3 Focal Plane Crystal Spectrometer 


The hrst high resolution observations of cooling clusters comes from the Focal 
Plane Grystal Spectrometer (FPGS) on the Einstein observatory. This instru¬ 
ment used a crystal to disperse X-rays according to the Bragg condition. It 
therefore could produce a high resolution spectrum by scanning in angle and 
therefore scanning in wavelength. This worked for an extended source like a 
cluster as well. 


Unfortunately, there were only a limited number of these observations and 
the count rates were not particularly hig h. These observations demon strated 
the existence of line emission from 8 ions ( Ganizares et ah 19791 1982[l in four 
different clusters. Some of these observations are contradicted by Reflection 
Grating Spectrometer observations as will be described in the next section, 
so it is likely that these observations were partially compromised by high 
background rates and, in some cases, relatively few photon counts. 


5-4 Reflection Grating Spectrometer Observations 


Observations made at high spectral resolution with the Reflection Grating 
Spectrometers on XMM-Newton have greatly clarihed the observational inter¬ 
pretation of the soft X-ray spectra of cooling clusters. Although some informa¬ 
tion can be gained from the exact shape of the Fe L complex with low spectral 
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Fig. 6. The cooling time as a function of radius for a sample of relaxed clusters as 
measured with Chandra. If left undisturbed, all pl asma within 100 to 200 kpc would 
have sufficient time to cool. Figure adapted from Voigt fc FabianI (2004). 


resolution spectrometers, RGS observations resolve emission lines from each 
Fe L charge state (Fe XXIV through Fe XVII) and therefore can place crucial 
constraints of the amount of gas present between temperatures of 0.4 and 4 
keV. 


The initial application of the standard c ooling flow model to a high resolution 


X-ray spectrum is shown in Figure 8 (jPeterson et al.M2nni[l . The top panel 


shows the spectral prediction between 5 and 35 A. Emission lines are expected 
from Fe L ions in roughly equal strength between 10 and 18. Particularly 
prominent are emission lines from Fe XVII and 15 and 17 Angstroms. These 
emission lines are produced primarily from plasma between 300 and 700 eV 
and represent the last major emitting ion before a cooling plasma would cool 
to sub-X-ray temperatures. 


High resolution spectroscopic evi dence for lack of coo l gas has been docu¬ 
mented in a number of clusters (IPeterson et ah 2001: Tamura et ah 2001a : 
Ka.astra et IDboOltlS akelliou et al •1120021: 1 Tamura et alJboOlb ll. Characteristi- 
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Fig. 7. The average temperature in radial bins for a sample of relaxed clusters 
of galaxies as measured with Chandra. The temperature and radius are scaled to 
^’ 2500 - A typical cluster shows a clear decline in the average temperature at the 
center, which is in agreement with many spectroscopic studies over the previous 
two decades. Figure is adapted from ( Allen et al. 2nnicli . 


cally, gas appears to be missing at near a third of the maximum te mperature 


This p henomena has also been documented in elliptical galaxies. IXn et al 


f)2nn2ll found that Fe XVII and Fe XVIII was present in the nearby ellipti¬ 
cal galaxy, NGC 4636. O VII, however, was not detected and has not been 
detected in any galaxy, group, or cluster core. 


A more detaile d anal ysis of the nature of the cooling flow problem is found in 
Peterson et ahl ()2nn3n . A sample of 14 clusters were analyzed in a uniform way 
to demonstrate that the cooling flow problem manifests itself at a fraction of 
the maximum temperature in the center. In this sample it is shown that hot 
clusters (5-10 keV) generally only show evidence of Fe XXIV-XXII emission 
and no other Fe L charge states. Intermediate temperature clusters (2-5 keV) 
show evidence for Fe XXIV-Fe XIX emission, but no Fe XVII and Fe XVIII. 
The coolest clusters and groups (less than 2 keV) show evidence of the entire Fe 
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L series but anomalously low levels of Fe XVII emission when compared with 
the standard cooling flow model. All of these observations are in contradiction 
with the standard cooling flow model, and suggest qualitatively that the model 
is violated at characteristically a fraction of the temperature and is more 
inconsistent with the model at lower temperatures. 


In this quantitative analysis ( Peterson et ah 200311 . the cooling flow model was 
separated into four different temperature ranges as shown in Figure 9. Then 
the normalization of each temperature range was adjusted to empirically fit 
the spectrum shown in Figures 10, 11, and 12. The results from this study are 
shown in Figures 13 and 14. These graphs plot each of the four luminosity de¬ 
tections for each clusters. The lines connect the points that are actual spectral 
detections and not just upper limits. The standard cooling flow model would 
be if all the points would lie on a horizontal line. A comparison of Figure 13 
and 14 suggests that when plotting the differential luminosity vs. the fraction 
of the temperature as in Figure 14, that the trend is more systematic. There¬ 
fore, the emission of a cooling flow is consistent with a power law in differential 
luminosity with exponent between 1 and 2. The cooling flow model, however, 
is strongly ruled out in many systems. 


At this point, there has been no significant evidence for any significant differ¬ 
ence between cooling flows of similar cluster mass. This poses some challenges 
for heating models where there is a time-dependence to the process. There 
also has been no evidence from X-ray spectroscopy that there are signihcant 
quantities of plasma at low X-ray temperatures above the expectation of the 
cooling flow model. In other words. Figure 14 shows a monotonically decreas¬ 
ing violation of the cooling flow model and no evidence for gas piling up at 
some intermediate temperature. This presents some challenges for continuous 
heating model that do not reheat the gas completely. For both of these rea¬ 
sons, it is not straight-forward to interpret these results in the context of any 
heating model regardless of the mechanism. 


Finally, it is worth noting that these results generally are displaying the X- 
ray spectrum of the entire cooling flow volume. There is considerable work 
that is continuing to actually study the changes in the emission spectrum as a 
function of spatial position inside the cooling flow. These studies are difficult 
since determining a spatially-resolve differential emission measure require both 
high angular resolution and high spectral resolution. Following, we discuss 
some attempts to study the spatial-dependence of the X-ray emission of cooling 
flow at moderate spectral resolution. 
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5.5 Recent Spatially-resolved Spectro-photometric Observations 


The RGS results obtained at high spectral resolution have been augmented 
with further studies at low resolution that attempt to study the spatial distri- 


bution in greater detail (iDavid et al.l2nni: Bdhringer et al.l2nn 

;Mo] 

endi & Pizzolatol 

2001: Schmidt et ah 2001: Ettori et ahl 20021: l.lohnstone et al. 

2001 

). In gen- 


eral, it has also been demonstrated that the low spectral resolution observa¬ 
tions obtained with XMM-EPIC and Chandra, which are far more numerous, 
are generally consistent with the RGS observations. Signihcantly smaller mass 
deposition rates have been measured when the cooling flow model has been 
applied. Similarly, cut-off cooling flow models have been shown to be consis¬ 
tent with the data. Needless to say, a larger range of models can be applied to 
the low resolution data and still be statistically consistent with the data. In 
general, it appears that there is no strong evidence for any significant amount 
of cold gas (below | of the maximum temperature) in any cluster. 


The radial prohle of the temperature distribution in cooling flows is currently 
in dispute. It is clear that the ICM has a much narrower range of temperatures 
than compared to the standard cooling flow model as established by the RGS 
observations. It is not so clear, howeve r , how narrow that distribution is at 
any given radius. Molendi fc Pizzolatol (2001) have argued that the ICM is 
nearly isothermal at a gi ven radius and any depa rture from this i s due solely 


to azimuthal variations. Buote et ah (200^ and Kaastra et ah (2004) have 


argued, however, that nearly a factor of two range of gas temperatures exists 
at any radius. 


X-ray cavities have been discovered in a number of X-ray imaging observations. 
The cavities are presumabl y excavated by cosmic r ays produced in outbursts 
by the central radio source (jMcNamara et al.ll200lh . These cavities are impor¬ 
tant in the study of cooling flows for a number of reasons. First, they represent 
direct empirical evidence that the cores of cluster have been disturbed. The 
work done by expansion of the bubbles assuming all the energy is deposited 


in th e cooling flow is within a factor a few to the value required ([Birzan et ah 


200411 . Second, they are relatively intact and coherent structures indicating 


the level of suppression of thermal transport processes possibly by a magnetic 
held. Therefore, they provide a source of external energy to the IGM and de¬ 
tailed study will eventually give us a complete picture of the transport of that 
energy. 


Gonstant pressure surface br ightness discontin uities, or “cold fronts”, are present 
in many clusters of galaxies ( Markevitch 2000fl . Gold fronts show that there are 
vast quantities of plasma oscillating with respect to the gravitational poten¬ 
tial. Some may be cooler, denser groups that have fallen to the center, whereas 
others may be cooled gas that already was present in the core of the cluster 
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but i s oscillating because of a merger in the cluster outskirts ([Churazov et al. 
200311 . 


Fine r ipples have been identihed in the cluster gas of Perseus fjFabian et ah 


2003all . These features are extremely weak and more study is needed to mea¬ 


sure the temperatures across the discontinuities and confirm the ir properties. 


Dete ctions of thes e features have also been claimed in Hydra A fjNulsen et al, 
2004 1 and M 87 ( Forman et al. 2003ll . It is unclear if these are an entirely 


separate phenomena from cold fronts. Given their proximity to the center, 
Fabian et al. (l2003a l suggested t hat these were weak shocks and sound waves. 
In addition. ISanders et ah ( 2004 1 have identihed an asymmetric “swirl” in the 
temperature distribution in the Perseus cluster. This observational feature is 
may be due to the effect of off-center cluster mergers. 


5.6 Observations of Cataclysmic Variables 


Although not the subject of this review, it is certainly worth noting that the 
X-ray spectrum of some binary accreting white dwarf systems, cataclysmic 
variabl es, actually resemb les the X-ray spectrum of the standard cooling how 
model ( Mukai et al. 2003ll . The how behind the standing shock in an accret¬ 
ing white dwarf involves the same physical processes as those expected for a 
cooling how in a cluster although the geometry may be diherent. Although 
there are a great number of diherences between clusters of galaxies and cat¬ 
aclysmic variables, these observations can be used to conhdently say there 
is no mistake in applying a model for a radiative-cooling dominated plasma. 
Furthermore, this points to no major difficulties in the atomic physics models 
or subtle problem with the plasma physics arguments we have used. 


5 .7 Definition of the Cooling Flow Problem 


We now briehy discuss what we believe the cooling how problem is, and how 
it might be resolved. Clearly, the proble m is quite complex and it is difficult 
not to see the problem in historical terms. Peterson et al. ( 200311 . for example, 
discussed a diherence between the soft X-ray cooling-flow problem and the 
mass sink cooling-flow problem. The former refers to the recent discrepancy 
seen in the soft X-ray spectrum between what was predicted and what was 
observed. The latter refers to the difficulty in detecting any by-products in 
cooling clusters from the hypothesized cooling-how plasma. 


These dehnitions, however, might just categorize our ignorance of the solution 
to the problem. The major difficulty is that: 1) the cluster plasma loses energy 
by emitting the very X-rays we detect, 2) efficient and distributed heat sources 
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are difficult to construct, 3) the cluster plasma appears to cool most of the 
way, but 4) evidence for complete cooling is utterly lacking. The cooling-flow 
problem as we see it is to understand what happens in the middle of that 
process. After examining whether cooling flows are ruled out, we discuss many 
ideas that might alleviate the cooling-flow problem. 

Note, from Fig. 6, that time variability on intervals longer than 10® yr cannot 
be the solution. 
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Isoboric Multiphase Cooling Flow Model 



' X ' 

■ 

2.5 

— 1 

kT^= 8 to 0 keV “ 


- 5 5 

1/3 Solor Abandonees I 

2.0 

z 4» X X 

— >- - 

_ 


— w A O s;- 4) 

^ XX ^ Li. ^ Li. 

X XX f X X 

- 

1.5 

^ 0) 0) . 0) 0) 

_ S LI.L. L. L. 

0 VIII J 

1.0 


-rryT 

lA N 

lA 0 

0.5 

I_ Bremsstrohlung Continuum 

.__ 


" 

. 


10 


15 _ 20 

Wovelength (A) 


25 


30 


Abell 1835 ond 2300 Mg yr ' Cooling Flow + kTj=8 keV Ambient Component 



15 20 

Rest Wovelength (A) 

Abell 1835 ond Empiriool Model with kTj=2.7 keV Cut Off 



15 20 

Rest Wovelength (A) 

Fig. 8. Three panels showing the cooling-flow model applied to the X-ray spectrum 
of the putative massive cooling-flow, Abell 1835. The top panel shows the standard 
isobaric cooling flow model (see the previous figure). The middle panel shows the 
model (blue) and the data obtained from the Reflection Grating Spectrometer on 
XMM-Newton (red). The model is clearly inconsistent with the X-ray spectrum, 
particularly in the prediction of Fe XVII emission line blends at 15 and 17 A. The 
bottom panel shows the cooling flow model compared with the data, except all 
emission coming from temperatures below 2.7 keV is suppr essed. The explanatio n 


for the success of this model is not known. Adapted from Petersou et al. I (l2nnih . 


The spectrum is taken from a 5 by 20 arcminute region of the core. 
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Isobaric Multiphase Cooling-Flow Model 



Fig. 9. The isobaric cooling flow model divided into the contributions from various 
temperature ranges, 0.375 to 0.75 keV (blue), 0.75 to 1.5 keV (green), 1.5 to 3 keV 
(yellow), 3 to 6 keV (red). 
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Fig. 10. Two models compared with the Reflection Grating Spectrometer data of a 
sub-sample of massive cooling flows. The blue histogram is the RGS data, the green 
curve is the standard cooling-flow model, and the red is an empirical model where 
emission from the standard cooling-flow model is allowed to be adjusted in specific 
temperature ranges. The spectrum is taken from a 5 by 20 arcminute region of the 
core. 
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Fig. 11. The same as the previous figure, but for a sample of intermediate mass 
clusters. 
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Fig. 12. The same as the previous two figures, but a sample of low mass clusters 
and groups of galaxies. 


42 


10"^ Photons/cmVs/A 









































s 

»f5IN 




HH 


Cooling Luminosity Limits 
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Fig. 13. The relative amount of emission in various temperature ranges from 
a sample of 14 clusters obtained with the Reflection Grating Spectrometer on 
XMM-Newton. Each cluster is a different color, and each cluster has four points. 
Points that are not upper limits are connected by a straight line. The standard cool¬ 
ing flow model predicts the same amount of emission in each temperature range or 
a horizontal line. The data, however, are clearly inconsistent with that model. The 
detected emission has a much steeper distribution in temperature a nd many upper 


limits are a factor of several below the prediction. Adapted from iPeterson et al 

( 2nn3h . 
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Cooling Luminosity Limits 



Fig. 14. Same as the previous figure, except the points are normalized to the max¬ 
imum temperature in the fit. A more systematic trend is seen where the emission 
follows a linear or quadratic form in temperature (dotted lines). There is significant 
scatter in the points, but there does appear to be a “self -similar” violat i on to the 
standard cooling flow model (dashed line). Adapted from IPeterson et~ (j2nn.'lll . 
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Fig. 15. High resolution X-ray image of the Perseus cluster of galaxies. The bright 
source in the center is the active galactic nucleus of NGC 1275. The two adjacent 
holes in the X-ray emission as well as the mushroom shaped depression in the upper 
right are believed to be cavities in the ICM that have been excavat ed by cosmic 


rays ex pelled from the active galactic nucleus. Figure is adapted from iFabian et al 
( 2nn5bli . 
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Fig. 16. Image of the Centaurus cluster of galaxies ov erlayed with Ha emis sion 


showing a filamentary structure. Figure is adapted from ICrawford et al 
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6 Are cooling flows ruled out? 


In summary, problems with the simple cooling flow picture emerged from two 
fronts. The first is the enormous implied mass deposi tion in some obiects where 
for example M ~ 1000 Mq yr“^ (e.g. PKS 0745-191: Ih^bian et a~ 1985li which 
is then not seen in cooled form (the ’mass sink’ problem). The central galaxy 
of this, and many other brightest cluster galaxies where the surrounding gas 
have tcooi < 3 Gyr, has much exc ess blue light, optical/UV/IR emission-line 


nebulosity ( Crawford et ah 19991 and refs, therein) and even molecular gas 


( Edge et ahlbooi h The total star formation rate and gas mass is often one to 
two orders of magnitude below such a large value of M operating for several 
Gyr. The flow could be somewhat intermittent but not overly so or short 
central cooling times would not be so com mon dPeres et ah jT98 h A non¬ 
standard I MF for the stars is one possibility ( Fabian et ah I ll982ll . but has to 


be extreme (jPrestwich et al.lll997li . 


The second problem is that a sim ple cooling flow spectrum is a poor fit to 
the data in the so ft X-ray band (jwhite et ah 1991 : Johnstone et ah 19921: 


Fabian et al.lll984li . Wha t was usually done wa s to fit the model spectrum of 
an isobaric cooling flow ( .Tohnstone et ah 1992[i plus an isothermal spectrum 
to represent the luminosity from gravitational work done, to the data from 
the inner regions. A deficit in soft X-ray emission was found in the data when 
compared with the model. This was then modelled as photoelectric absorption 
due to an absorber intrinsic to the central regions. The lack of any obvious 
absorption in very soft ROSA T spectra was attribu ted to emission from the 
intervening gas in the cluster ( Allen fc FabianI 1997 ). 


Photoelectric absorption i ntrinsic to the cluster d oes not however appear in the 


(Chnrazov et al. 20081. Nor is it api 

parent in the detailed XMM/RGS spectra 

of cool cluster cores (I 

Peterson et al. 

2001 :J 

^amura et al.ll2001al: Kaastra et al. 

2001 : Sakelliou et al. 

2 OO 2 I: IPeterson et al. 

20081. Absorption has generally 


been abandoned as an explanation, although no detailed test has been made 
of a model in which the absorption is intimately linked with the coolest cooling 
gas clouds. Note though that the presence of warm emission-line nebulosities. 


dust- lanes and clear optical absorption lines (jCarter et al.lll997t ISnarks et af 


199711 show that at some level there must be distributed intrinsic absorption. 


Curiously, the ’missing’ soft X-ray luminosity (i.e. the above soft X-ray deficit) 
is close to the luminosit y emitted from the emission-line nebu losity commonly 
found in these systems ( Fabian et al 1l2nn2al: ISoker et al 1l2nn4 . This has led to 
the suggestion that the gas is cooling but not radiatively once its temperature 
has dropped to ~ 1 keV or so. It c an be cooled by mixing ( Fabian et al 1l2nn2ah 
or conduction ()Soker et al.N2nn4ll with cold gas. A possible picture would be 
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that the gas cools from 100 to within 20 kpc radiatively as a single-phase 
flow by which radius its temperature has dropped by a few. It then shares its 
temperature with cooler gas and thereby rapidly drops down to below 10® K. 
The cooler gas is close to the peak of the cosmic cooling curve and radiates 
the energy in the UV and optical bands (and IR if there is molecular gas and 
dust). Such a solution to the soft X-ray dehcit still suffers from the mass sink 
problem. 


A further situation which could give a soft X-ray dehcit is if the metals in the 
ICM are highly inhomogeneous. They are presumably introduced in a very 
localized manner and at high abundance by stars and supernovae. If they do 
not mix but after time just reside at a similar temperature to their surround¬ 
ings then when bremsstrahlung cooling dominates (which it does above about 
1 keV for Solar abundances) then all the gas cools together. But when it has 
cooled so that line cooling dominates, the high l y enriched clumps woul d then 
cool rapidly and drop out (IPabian et ah 2001: Morris fc Fabian 200311 . This 
has been considered in part bv iBohringer et ah ( 2002il . although no conclusive 
test has been performed. 


The main reason that complete cooling flows, in which gas cools from the 
cluster virial temperature to well below 10® K, are often now considered to be 
ruled out is the high spectral resolution XMM/RGS data ( Peterson^^^^ II 2 OOII 


Tamnra et ah 2001a : Kaastra et al. boOltl^kellion et al. ll2002l:IPeterson et af 


20^ h Chandra spectra lead to a similar conclusion ( David et al •l200ll: l Allen et ahl 


2001all . These spectra have been covered in depth earlier in this review. We 


summarize here to note that most data are consistent with the ICM temper¬ 
ature dropping within the cooling radius by a f actor of about 3 (some times 
signihcantl y larger values are f ound e.g. A2597 ( Morris fc Fabianll2005ll and 
Centaurus ( Fabian et al.ll2005a l. Data from all wavelengths are consistent with 
a residual flow ranging from a few Mq yr ~^ to atypical value of tens M^t^ yr ~^ 


up to ~ 100 Mrr^ yr ^ in some cases 


(e.g 


Baver-Kim et al. ( 2002il : Wise et al 


( 2004 1: McNamara et al. ( 2004ll i. These values are generally compatible with 


the ongoing star formation seen. 


6.1 Cooled gas and star formation 


The brightest galaxy in X-ray peaked clusters often has excess blue light indi¬ 


cating massive star formation and IR/optical/UV emis sion lines (jCrawford et al. 


I 999 I: Donahue et al. 2000 : Hicks fc Mnshotzkv 2005ll . The regions are often 


dusty so the determination of the star formation rate depends on uncertin 
dust corrections. It can be tens Mq yr“^ to about 100 Mq yr“^ in the case 
of A1835. In a typical BCG it is one to ten per cent of the mass cooling 
rate inferred from the hotter X-ray emitting gas assuming it cools completely. 
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UV-excess BCGs probably have the highest rate of star formation of early- 
type galaxies at the present epoch. The brightest galaxy in non-X-ray-peaked 
clusters do not have this emission. 


The excitation of the emission-li ne nebulosity, which o ften has a hlamentary 
structure (e.g. the Pe rseus cluster: Conselice,^^^ Ibnnnl: the Centaurus cluster: 
Crawford et ah 2005; the Virgo cluster: Sparks et ah 12^004 1 is likely related to 
the star formation but in detail there are problems. Hotter emission is required 
than expected from the stars observed and there are no stars obvious within 
the outer hlaments. The gas, even in the outer hlaments , contains dust and 
molecular hydrogen ( Hatch et al 2005a : Jaffe et ah 2005ll . Simple models for 
the molecular gas imply higher p ressures than for the surrounding ga s and 
very short radiative cooling times ( Jaffe et al. 2001 : Wilman et al. 2002ll . The 
hlaments may have been pulled out from a central reservoir of cold and warm 
gas by the action of radio bubbles ( Fabian et al. 2003bl: Hatch et al 2005bll . 
Where the central gas comes from is not clear although radiative cooling is 
likely. A residual cooling how is therefore taking place in many clusters. 


Gas at intermedi ate temperatures of a bout 10^’® K is seen through OVI emis¬ 
sion with FUSE. Oegerle et al. ( 200 found suc h emission from A259 7 and 


recently it has been found in A426 and A1795 ( Breeman et al. 2005ll . The 


inferred cooling rates in the 30 arcsec FUSE aperture are 30 — 50 Mq yr 


-1 
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7 Heating 


Some heating is alwavs expected in the central regions of clusters. Examples 

are super novae (IS 

ilk et ah Il986l: Domainko et ah 20041). an active central nu- 

cleus (lBailevlll982 

: Tucker & Rosneilll983HPedlar et ah 1990l Tabor & Binnev 

lOOOl: Binnev & ’ 

aboilI 1995 I) and manv more recent papers cited in Sec 7.2- 

7.5), conduction ( 

lakahara & Takaharal19791 Binnev & Gowielll981: Stewart et ah 

1984: Friaca 1986 

Bertschinger & Meiksin 1 19861 iRosner & Tuckei 19891) and 


many more recent papers cited in Sec 7.1). 


A problem with heating the gas is that the cooling rate is proportional to the 
density squared whereas most heating processes are proportional to volume. 
This tends to make the gas unstable and means generally that the cooler 
denser gas will carry on cooling while hotter surrounding gas heats up. The gas 
appears to cool by about a factor of three and then stop cooling. A mechanism 
to do that is not obvious, since the gas does not appear to be piling up at 
the lower temperature. Indeed it seems that the gas temp erature prohle is 


’frozen” and has been so for some Gyrs ( Bauer et ah 20051) . 


The prohle of fcooi is similar in many clusters with a common central minimum 
value for fcooi of about 200 Myr (Fig. 6). This strongly suggests that heating is 
continuous, at least on timescales of 10^ yr or more and is spatially distributed. 
Moreover, no shock waves are apparent in these regions so any mechanical 
energy injection must be subsonic. 


Brighenti fc Mathews ( 20031) (see also Ruszkowski fc Begelman 2002|) provide 


several ID simulations of a cluster core in which various levels of heat are 
injected at various radii. The results which best ht the temperature prohle 
of a real cluster are those without heating. This demonstrates that a heating 
model is not entirely straight-forward. 


Some authors use a combination of processes. Although several are operating 
in any cluster, it is unlikely that they have similar weight, especially when it is 
recalled that the cooling how problem exists in objects ranging from elliptical 
galaxies, through groups to the most massive clusters. It applies to a range 
of about 100,000 in X-ray luminosity and over 15 in temperature. We are 
therefore seeking a wide ranging, quasi-continuous, gentle, distributed heating 
process 


We now consider conduction and AGN heating in detail before discussing other 
models. 
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7.1 Heat conduction 


The outer cluster atmosphere beyond the cooling radius represents a vast 
reservoir of thermal energy. Conduction can transfer heat into the cooler cen¬ 
tral regions. The rate of conduction is however highly uncertain due to the 
presence of magnetic fields which are probably tangled. 


Most workers have adopted Spitzer ( Soitzer et ah 1962il conductivity (ks cc 
which is appro priate for an unmagn e tized gas, and then assumed a 


suppression factor /. iNaravan fc Medvedev! ()2nnill suggest that / is about 
1/5 since held line s wander away from each other exponentially when tangled 
fsee iTribble 19891 for another discussion). 


Comparison with the data initially looked promising (jZakams ka fc N aravanI 
29981:IVoigt et al.l2nn2h then u nder further inspec tion ruled it out ]k im fc NaravanI 
2991 ia yoigbjzjabi^ 2994). Voigt et ah ( 2992fl (see also lKaastra et ^l2994l 
and Voigt fc Fabianll2994ll used the density and temperature prohles to deduce 
the heat transfer rate within the gas and thus the effective conductivity 
required. For several high temperature clusters, Keg < ks , but for clusters 
where most of the gas is below 5 keV then > d^s. T hen seems unlikely, 
however, but this may be possible ( Cho fc LazarianI 2994 1. 


Cho et ahl (j2998ll proposed that turbulent conduction might work, whereby 
large subsonic eddies cause heat to be transferred i n radius. This does ap¬ 
pear to operate fast enough to explain the data ( Kim fc NaravanI 2993a : 
Voigt fc FabianI 1299411 provided the gas is highly turbulent. We shall return 
to the issue of turbulence later when discussing heating by the central radio 
7.2). 


source 


We note that conduct ion as a solution has a fine tuning problem fjNnlsen et al, 
19821: Bregman 1988ll . If the gas starts out isothermal, then either the gas 


cools in which conduction cannot stop it (conduction has a steep positive 
temperature dependence) or conduction does work and there is no cool gas. 
How a prohle resembling a cooling flow would occur is difficult to understand. 


7.2 Heating by a central radio source 


Most central cluster galaxies where in the surr ounding gas is less than a 
few Gyr have a nucleus radio source (|Bnrng 1999ll . The radio luminosity of 
these sources has a wide range from the powerful FRII objects Cygnus A and 
3C295 to weaker FRI objects such as M87. ROSAT HRI data clearly showed 
that the FRI lobes of 3C84, the central radio source in the Perseus cluster. 


have displaced the ICM and so are strongly interacting there (jBohringer et al 
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Fig. 17. The coefficient of electron thermal conduction that would be required to 
balance X-ray cooling in the cores of a sample of clusters of galaxies. These data 
are derived by measuring a mean temperature and density as a function of radius 
and calculating the conduction coefficient to cancel radiative cooling. Each color is 
a different cluster and each point is a measurement at variou s radii. The lines ar e 


prediction for the magnetic field free conduction coefficient of iSnitzer et al.l (|1 


It is unclear what to expect for the effective conductivity in a cluster of galaxies with 
a complex magnetic held topology (see text), but it would probably be below the 
Spitzer value. These results probably indicate the conduction is not strong enough 
to cancel cooling by itself, but the value required curiously happens to be within an 
order of magnitude of the canonical value. 


I 993 I 1 . The hi gher spatial resolut i on of Chandr a has revealed further details 


of that source ( Fabian et ah 20091 2003a, 2n05bll and enabled the discovery of 
many more, similar, interactions. 
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Fig. 18. The total heating power required to balance radiative cooling a sample of 
cooling-flow clusters over a Hubble time. The left hand side also indicates the mass 
of the black hole required if ten perce nt of its rest mass w ere applied directly into 


heating the ICM. Figure adapted from Fabian et al. ( 20021 1. 


7.2.1 Overview of radio bubbles 


Holes in the X-ray surface brightness coincident with radio lobes are commonly 
seen and generally referred to as radio bubbles. They are presu med to be rela¬ 


tively empty of normal thermal gas (see limits by Schmidt et al lSchmidt et al. 


200111 and mostly hlled with a plasma consisting of relativistic electrons, pro- 


are sc 

Centa 

A259' 

en in Hydra A (McNamara et ai.l 2000ll: A2052 ( 

31anton et al.l 

2001 

); 

1 ; 

urns (IPabian et al.ll2000ll: M87 (lYonne; et al. 20021: 

Forman et al. 

2003 

7 (McNamara et al.l200ll: Clarke et al. 11200511 and many more (IBirzan et al. 

2004: Dunn & Fabian 2004 Dunn et al.|2005|l. 


Sometimes outer bubbles are seen (the NW on e in the Perseus cluster was hrst _ 

seen in images from the Einstein Observatory, ( Fabian et al IhOSlI: iBrandnardi-R avmont et Id. 


198111 which are considered to be ghost lobes. Synchrotron and other losses 


have depleted the population of radio-emitting electrons to the extent that 
they are not detectable in the radio (particularly not at high frequencies) yet 
the other electrons, protons and magnetic held in the cavities still exclude the 
X-ray emitting thermal gas. Weak, low-frequency radio emission pointing at 
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Fig. 19. The instantaneous heating rate required as a function of radius for a sam¬ 
ple of massive cooling-flow clusters. These results indicate that the heating has to 
be spatially dis t ribute d over a substantial physical volume. Figure adapted from 
Voigt fc FabianI ( 20041! . 


the o uter bubbles in the Perseus cluster supports this hypothesis fjFabian et ah 


20021! 


The bubbles discussed so far have radii from about 1 to 15 kpc. Recently, giant 
bubbles of radius 100 kpc have been found 100-200 kpc from t he centers of the 
Hydr a-A nucleus ( Nulsen et al ]l2004l! and in MS07 35.6-1-7421 (i McNam ara et al 


2005l! . and some bubbles with little energy (e.g. Blanton et al. ( 2004l! !. The 


minimum energy of a bubble can be estimated from its volume V and sur¬ 
rounding pressure P as T^bubbie = PV■ If the interior pressure is entirely rel¬ 
ativistic in gas with a ratio of specihc heats of 4/3 then T ^ K„bhiP = 4P1/. The 


energies are signihcant for typical bubbles (jBirzan et al 
106 i gj-g fQj. giant bubbles. 


120041! and exceed 


7.2.2 The simple theory of bubbles 


Bubb les from radio sources were hrst predicted by Gull & Northover ()Gu 


973l! and st udied ana lytically before th e Ghandra/XMM- Newton era b y 


l&^orthoveij 


( 2000 , 200ll! (see also Soker et al. 20021! and simulated by iHeinz et~ah (1 


Ghnrazoy et al 


They are inflated by jets from the nucleus and typically form in pairs either 
side of the nucleus and seemingly attached to the nucleus. If a jet of power L 
steadily inflates a bubble of radius R then, after an initial supersonic phase 
and assuming radiatiye losses are negligible, the bubble expansion with time 
t follows from the conseryation of energy as Lt = PV oc PR^. Thus R oc tl/3 
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and the expansion speed of the bnbble ^bubble oc The bnbble is of mnch 

lower density than its snrronndings so is bnoyant and rises at a velocity close 
to the local gravitational free-fall valne Wgrav When ^bubble drops below Vgrav 
the bnbble detaches from the jet and rises. A new bnbble then forms if the jet 
remains powered. 


A simple estimate of the power injected into the bnbbles can be obtained by 
dividing T^bubbie by an estimate of the age of the bnbble. which can be obtained 
from the above sc aling (say, R/v^r&v) (iBirzan et al.l2004 lTA inn k Fabianll2004 
Dnnn et al. 200511 . These are typically several orders of magnitnde higher than 


the inferred radio Inminosity of the jets and nnclens showing that the radio 
Inminosity is a poor gnide to the power of a jet. Snch radio jets and lobes 
can have a radiative efficiency (ratio of observed synchrotron Inminosity to 
mechanical power) as low as 10“^ or even 10“^. This, together with time vari¬ 
ability of the jets, indicates that the very poor correlation b etween the radio 
power of the central sonrc e in a clnster and the heating power ([Voigt fc Fabian 


2004 iKaastra et al.ll2nn4ll does not necessarily imply that the mean jet power 


does not correlate well. 


The power of the jets is mostly snfScient to offset r adiative cooling in the 
clnster core provided it is steady. iBirzan et ahl ([20041 did however note that 
it is insnfficent in some objects. This conld indicate that radio sonrce heating 
is not a general process or that it can vary by factors of a few on the bnbbling 
timescale, which ranges between ~ 5 Myr and ~ 50 Myr. This timescale 
is still shorter than the inner radiative cooling timescale which is typically 
100-300 Myr. 


A farther problem they raise is that bnbbles were only fonnd in 10 per cent of 
the clnsters they examined, so they either do not occnr in many clnsters or if 
they do then they do not occnr often. This resnlt could however be a selection 
effect from the examination of a sample of clusters mixed in type (i.e. whether 
the cooling flow problem applies to that cluster or not) and signal-to-noise. 
A simple inspection by ACF of Chand ra images o f the 3 0 per cent of the 55 
brightest clusters in the Sky listed by Edge et al. (199^ which have central 
Aooi < 5 Gyr shows that all but one have clear bubbles, so we suspect that 
bubbles are indeed common enough to be a viable ingredient in th e solution to 


the co oling flow problem. A more detailed analysis of the issue bv iDnnn et af 


( 2nn5ll shows that at least 75 per cent of the cluster cores needing heating 
have bubbles. This means that the duty cycle of bubbles is such that they are 
(detectably) present for 75 per cent of the time. 


Finally, before examining numerical simulations and then the heating mech¬ 
anism in detail, we note that the total energy required to stave of cooling 
for several Gyr in a luminous cluster implies a large central black hole mass 
( Fabian et al. 20021: assuming a mass-to-energy conversion effiency of 10 per 
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cent). iFuiita fc Reiprich ( 2004) compared bl ack hole masses predicted by the 
velocity dispersion relation ( Tremain3 20021) and found that they often fall 
short of what is required. However for the very few systems where the mass 
has been measured (e.g. M87) there is no problem. Perhaps the conclusion 
here would be that massive black holes in cluster cores lie above the Mbh — <7 
relation for lower mass galaxies. 


1.2.3 Numerical simulations of radio bubbles 


Many groups have now carried out nun ierical hydrodynarnical s i mulations of 


the 

)ehaviour of bubbles in the ICM (iBriiggen & Kaisei 2001: 

Ouilis et al. 

2001 

Briiggen & Kaiser 

2 OO 2 I: BruggeiJ2003l:lBasson & Alexander 120031: lOmma & Binnev 

2004 

Omma et al.l2004 

R.uszkowski & BegelmarJ2002:lR,uszkowski et al.ll2004allb: 

Robinson et al. 20041: Vecchia et al.l 2004: |Reynolds| 2005|). 


Most have been 2D or 3D and use the FLASH or ZEUS codes tsee lGardini fc Ricker 


20041) for a brief discussion and comparison of many of the simulations). Sev¬ 
eral assume that the surrounding gas is isothermal. Many produce bubbles 
which are very unstable (to Rayleigh-Taylor and Kelvin Helmholtz instabili¬ 
ties depending on the motion of the bubble at the start of a run) and collapse 
once they have risen more than their own height, which is most unlike the ob¬ 
served bubbles. Indeed most simulated bubbles do not look like the observed 


ones. 


Several of the early simulations relied on the dragging out of cooler gas from 
the smallest cluster radii by the bubbles as the explanation of why so little 
cool gas is seen in real clusters. This explanation however only works over 
the lifetime of the cluster (say 5 Gyr) if the cooler gas is dragged out beyond 
the cooling radius (i.e well beyond 100 kpc), otherwise it either falls back in 
or mixes in and reduces the cooling time of gas just beyond the center. A 
catastrophic cooling flow is only being postponed for a while. Although this 
explanation can work temporarily it cannot provide a comprehensive solution. 


When heating is estimated in simulations it is often unclear what it is due to. 
Some work is done by the buoyant bubbles rising and the sinking of cooler 
gas. This i s the heating n iechanism of the ’efferve scent heating’ approach of 
Begelman ( 20011) (see also iRuszkowski et al. l l2004^ . PdV work is done as the 


bubbles are made in the hrst place. That energy presumably propagates as a 
sound wave and need not be dissipated locally unless shocks are involved or 
the ICM is viscous. 


Heating by a mixture of weak shocks (at small radii ) and viscous dampin g of 
sound waves (at larger radii) has been proposed by Fabian et ah ( 2003bl) on 
the basis of the ripples seen in deep Chandra images of the Perseus cluster 
(see also Forman et al 2004 for a discussion of ripples in the Virgo cluster). 
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Viscosity has been i ncluded in simu lations bv iRnszkowski et alJ ()2nn4aybri and 
by Reynolds et al. ( Revnold'i 20051) and appears to dissipate 20-30 per cent 
of the PdV energy from the bubbles. Transport of the energy by sound waves 
and dissipating it over a lengthscale of 50-100 kpc provides a fairly gentle 
and distributed source of he at. The bubbles in a viscous medium also better 
resemble the observed ones ( Reynold^ 20051) . 


Fniita fc Suzuki ( 20051) and Mathews et al. ( 20051) argue that strong sound 


waves will shock and heat the innermost region and not deposit energy further 
out where it is need ed. However a very deep image of the Perseus cluster 
( Fabian et al. 2005bl) shows that the shock (Fig. 15) is isothermal. This may 
indicate that thermal conduction operates at least within the inner parts of 
the cool region to share the energy released by the bubbling process. 


7.3 Is the I CM turbulent? 


Several workers et ahllioO^ llnogamov fc Snnvae\( l2003l: IChandrai]l2004 l: 
Fniita fc ReinriclJ 2004 ) suggest that the ICM is turbulent. Such motions can 
transport heat and the dissipation of the turbulence is itself a source of heat. 


Schnecker et ahl ()2004l) have measured brightness and temp erature fluctua¬ 


tions i n the Coma cluster (which is not a cooling flow cluster). IVogt fc Ensslin 
()2004l) have determined a limited turbule nt-like spectrum in the Faraday Ro¬ 
tation Measures of the Hydra cluster core. Ensslin (1200, 5l) have extended their 
work to predict velocities and length-scales for the turbulence in a number of 
cool core clusters. They argue that a small-scale turbulent dynamo is main¬ 
taining magnetic helds and that turbulent dissipation can balance radiative 
cooling. 


Contrary indications have been identihed (Fabian et al 2003 b) in the highly 


exten ded, and often linear, optical hlaments in the Perseus (IConselice et al. 
2000 ) and other clusters (e.g A1795, Cowie et al ■119831: Cent aurus. lCrawford et al 


200,51). Such hla r nents have coherent v e locity helds and sm all velocity widths 
()Hu et al.lll985l: ICrawford et al.ll2005l : iHatch et a TI l2005bl) so must be rela¬ 
tively old (perhaps 100 Myr) and in thermal pressure equilibrium with the 
surrounding ICM. The optical surface brightness of such hlaments translates 
to a thickness of the emission regions of less than a pc. If the ICM is tur¬ 
bulent and pushing on a hlament, it will respond after interacting with its 
own column density which means a kpc or so of ICM, given the thousandfold 
higher density in the hlament. This is only a few arcsec in nearby clusters so 
the apparent linear nature of the hlaments would soon be destroyed by a fully 
turbulent ICM. They may however rehect streamlines and could have been 
drawn out by the motion of bubbles. Note the similarity of the horseshoe¬ 
shaped hlament in the Perseus cluster to that of streamlines behind rising 
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water b ubbles: see Fig. 184 in 

(l2nn3hli 


van 


Dvkelll982l and discussion bv iFabian et al.l 


The inner ICM in clusters with signihcant linear hlaments may therefore have 
subsonic random flows, particularly radial behind bubbles. Whether this is 
turbulence depends on how turbulence is dehned. It d oes not resembl e fully- 
developed hydrodynamical turbulence (see Fig. 186 in van Dvk3 1982ll but i 


IS 


not dissimilar frorn flow' s in some magnetohydrodynamical simulations 
Schekochihin et ah 2004 1. 


T-g- 


What is unclear is the driving mechanism for turbulence in a cluster core. 
A subcluster merger may well have made the Coma cluster ICM turbulent, 
but such mergers will n ot have strongly affected the dense core cluster under 
discussion here (but see lFuiita fc Reiori chi 1200411 . Driving turbulence into the 
dense core will require considerable energy. Rising radio bubbles are a possible 
driver but whether their motion leads to a turbulent cascade of energy to much 
smaller scales must be determined. The viscosity of the ICM may play a key 
role in this issue. 


The abundance gradients seen at the centers of many cool core cluste r s pro 


vides a further constraint on turbulence and diffusion. Rebusco et ah ( 2005ll 


have modelled the production and maintenance of the gradient in the Perseus 
cluster and find a diffusion coefficient D vi 


2 x10^® cm‘ s 


2 c-i 


Graham et ah 


( 209511 have carried out a similar analysis on the more abundance-peaked Cen- 
taurus cluster to obtain a value of D ~ 5 times smaller. 


7-4 Multiphase flows 


Large scale global flows have been considered tha t transfer signifi c ant therma. 1 
energy from one region of the cluster to another. Mathews et ah ( 20081 2004) 
have discussed flows which move in both radial directions. To be long lived, 
both mass and energy must be supplied to the inflowing gas over a large 
volume. The energy is assumed to be derived from bubbles. 


Fabian et al.l ()2003all considered a multiphase flow in which outer denser blobs 


fall inward and mix with inner, cooler blobs, thereby tapping the extended 
gravitational potential. Unless there is some outflow or cooling, however, mass 
builds up near the cluster center. 
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7 .5 Role of Magnetic Fields and Cosmic Rays 


A cooling flow will am plify tangled magnetic fields in the intracluster gas 
(jSoker fc Sarazinlll99fltl . Magnetic fields in the general ICM have been found 
throu gh radio observations of Faraday Rotation ( Kim et ai 1ll99l[lciarke et ^ 
200411 ■ Such observations also s how it rising in cool core clusters where the 


highest values have been found ( Carilli fc Tavloi 2002ll . The results indicate a 
roughly cell ular structure for t he magnetic held of a few kpc (e.g. the Centau- 
rus cluster, Tavlor et ah 2002ll . Such a struct ure may occur from a spectrum 
of magnetic helds perhaps due to turbulence (Ensslin 2005ll . Near the center 
the magnetic pressure can be 10 per cent of the thermal pressure, perhaps 
more. An important effect of the magnetic held is that it can dramatically 
reduce the microphysical transport processes in the ICM. 

Cosmic rays are a likely additional component in cluster cores. Relativistic 
electro ns produce the mi nihalo in the X-ray peaked region of the Perseus clus¬ 


ter (see ICitti et al.ll2004l for a theoretical discussion ) and may produce a hard 


X-ray component by inverse Compton scattering ( Sanders et ah 2005ll . The 


ehects of cosmic rays could be important in mediating the how of mechanical 
energy in the core. 


Cen (200^ notes that a cosmic ray phase can make heated gas thermally 
stable. A hoor temperature of 0.3 times the local ambient temperature is 
derived if the cosmic ray pressure is 1/3 of the ambient value. This situation 
is diherent from the similar factor seen in observations, which applies across 
the whole cool region, not just locally. Note that gas is not observed to ’’pile 
up” at any particular temperature. 


If cosmic rays are spread through the core in many small bubbles, then they 
would be very ehective in dissipating the energy in sound waves (Heinz & 
Churazov 2005). Such bubbles must be small not to be detectable in deep 
images (eg. Fabian et al 2005b). 


7 .6 Feedback 


If AGN heating balances cooling on timescales of ~ 10® to a few 10® yr then 
some feedback is needed to prevent either a cooling catastrophe (none has 
yet been found) or an event heating all the central gas and blowing it out. 
Bondi accretion onto a central black hole provides a link between the two 
regions but the radius range from the Bondi radius of say 50 pc to the ICM 
at 50 kpc is a factor of one thousand. The volume occupied by the black hole 
and its accretion radius which is required to provide the heating is one part 
in a billion! 
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Such numbers emphasize the problem face d by the necessary feedback. Some _ 

attern p ts to model feedback are discussed by Ciotti&_Ostrife J2OW Kaiser fc Binnev 
f 20031: Nulsen et al. (120041. Soker fc Pizzolato (12005 h lOmma fc BinnevI (12004 h 


Binnev pOOJl : Hoeft fc Briig-genl ( ~04h : ( Churazov et al. 12005 1: ( Pizzolato fc Sokei 

200,^ . The common occurrence of star formation in the central galaxy and also 


of giant bubbles shows that any feedba ck need not be perfect on timescales 
less than 10® yr. Donahue et al. ( 2005ll study two clusters which show little 
central temperature drop or central radio source yet have central cooling times 
of ~ 1 Gyr. Such clusters, which are in the minority, may be recovering from 
giant outbursts. 


7.7 Other heat sources and mechanisms 


Cosmic ray heatin g has been invoked in several models (jColafancesco et al 
2004 : Totanill2004 1. Generally to be effec tive, the pressure in cosmic rays would 
have to exceed that of the thermal gas ()Loewensteinlll99lll . Excess ionization 
by suprathermal electrons could cha nge the sp ectral appearance of gas in a 
cooling flow. This has been studied bv lOhl ( 2004^ . who hnds that such electrons 
provide r nuch more heating th an ionization, so reverting to the problems raised 
above bv iLoewensteinI ()l99lh . 


Hea ting by galaxy motions have been revived by |E 


also iFaltenbacher et al.ll2nn5ll following early work by 


-Zaute^^ f 2nn4ll 


see 


mmCm- It may be 


a viable distributed source of heat if the mass-to-ligh t ration of the galaxies 
passing through the inner regions exceeds 10 (but see iKimI 1200511 . 


Dark matter interactions have been discuss ed by several authors ()Qin fc Wu 


2001 : Totanill2004 : Ghuzhov fc Nuss ^l2004h . Similarities with the s olar 


corona 


have b een discussed bv iMakishimal ()2001rl and iKaastra et ahl (1200411 . iBriiggen fc Ruszkowski 
f 2OO5II discuss that viscous heating, if it occurs in real MHD plasmas, could 
provide signihcant heating during structure formation. 


One last possibility is that sedimentation a cts (IPabian fc Pringle 19771. He¬ 
lium then accumulates at the cluster center ( Gilfanov fc Sunvaevlll984 1 possi¬ 
bly explaining the puzzling low metallicity found at the center of some nearby 
clusters observed with Ghandra. The main low temperature cooling gas is then 
helium rich which can produce much weaker emission lines. Stars formed from 
helium rich ga s will have short lifetimes so avoiding some constraints on the 
cold mass sink ( Fabian et al 112003d: I lEttori fc Fabiar 11^ • Tangled magnetic 
helds will of course signihcantly reduce the prospect of any sedimentation. 
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8 Discussions 


The simplest explanation for the common appearance of cold core, X-ray 
peaked clnsters is that, when averaged over tens of Myr, the radiative cooling 
is balanced in part by distribnted heating. Thermal condnction as a means of 
distribnting heat from onter gas is rnled ont for low and intermediate tem- 
peratnre clnsters. It may however have a role in spreading the energy in the 
central parts. A plansible mechanism is the dissipation of energy propagating 
throngh the ICM from a central radio sonrce. Snch a process stems mas¬ 
sive cooling onto the BCG which wonld otherwise gain a total stellar mass 
S> 10^^ Mq. The process is th erefore a vital ingredient in stopping the growth 
of th e most massive galaxies ( Fabian et ah 2nn2bl: Benson et ahlboil^ Binnev 


200411 . NOte that most semi-analytic models for galaxy formation (e.g. Kanff- 


mann et al 1999) already needed to snppress cooling in massive haloes in order 
to match observation. 


Difhcnlties and donbts remain with regard to the issnes of the energy dissipa¬ 
tion and distribntion processes which are tied in with the transport processes 
in the gas. Similarly, it is not clear how the feedback manages to prodnce snch 
similar cooling time prohles in systems where temperatnres and thns masses 
differ by over an order of magnitnde? There still remains the possibility that 
some process not yet foreseen, or at least not well stndied will eventnally prove 
more important than the effect of the central radio sonrce, or will at least be 
important in mediating its effect. The central radio sonrce is so common and 
so energetic however that it mnst at least be part of the solntion. Similarly, 
the motions of galaxies or interacting dark matter, if it exists, conld be im¬ 
portant in heating clnster cores along with the AGN. Given the wide range of 
objects in which a balance is reqnired, we snspect that a single mechanism is 
dominant, rather than several. 


The need for a heating-cooling balance is in a time-averaged sense, over inter¬ 
vals of abont 10® yr. In most cases the heating has not been so energetic as to 
drive gas ont of the inner regions nor so weak as to allow mnch cooling at very 
high rates. Examples of objects at the extremes are Hydra A and Gygnus A 
where heating is high (bnt is dnmped mostly at large radii ontside the cool 
region) and A2597 and RXG1504.1-0248 where cooling appears to be high (in 
the latter obi ect over 70 per cent of the total X-ray Inminosity emerges from 
the cool core; Bohringer et al. 2005ll . 


In most objects residnal cooling at a rate of abont 10 per cent of the simple 
nnheated cooling rates appears to occnr. It conld be larger if non-radiative 
cooling, dne to mixing say, is occnrring. Stars form from the cooled gas giving 
the excess bine light seen in the BGGs. Mass loss from snch stars can make 
the cooled gas dnsty and radio babbles drag some of it back ont to large radii. 
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The time evolution of the heating / cooling balance is little understood. We 
suspect that the common temperature drop associated with the short central 
cooling times is due to radiative cooling and that heating only came into 
balance when the overall temperature structure was in place. Perhaps the 
central galaxy and its central BH grew until balance was achieved, and growth 
requir ed cooling. Compa rison of samples of clusters at z = 0.22 with ^ = 
0.056 ( Bauer et ah 2005ll shows little surprisingly change in the distribution 
of cooling times at 50 kpc. The results imply that any balance was established 
well beyond 2 ; rsj 0.3. 


Few massive cool core clusters are known at much higher redshifts than that 
of RXJ1347 at 2 ; = 0.44. This may however be a selection effect. They will 
be absent from X-ray cluster samples if the central BH is a bright X-ray 
source such as a quasar. If a central quasar outshines the host cluster in X- 
rays then the object will generally be classified as a quasar. H1821-1-643 is a 
good example of a bright quasy in a n X-ray peaked cluster at intermediate 
redshift {z = 0.297, Fang et ah 2002ll . Searches for clusters around powerful 
radio-loud quasars and galaxies have found some examples at 2 : = 0.5 — 1.1 


f Worrall et ah 2001 : Crawford fc Fabian 20031: Siemiginowska et ah 2005ll but 


no complete searches have been done. 
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9 Future work 


Further detailed deep studies with Chandra and XMM-Newton as well as 
future studies with Constellation-X and Xeus are needed in order to better 
understand the heating/cooling balance. There is considerable potential for 
studies of the cool and warm gas and dust around the BCG using HST, Spitzer 
and ground-based telescopes. Larger cluster samples are needed, particularly 
at medium to high redshifts. Determining the extent of a heating-cooling bal¬ 
ance (or not) in groups, elliptical galaxies, and galaxy bulges is also important 
as is the evolution of any balance in all massive objects. 
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